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We consider a classical random-bond Ising model (RBIM) with binary distribution of itK bonds 
on the square lattice at finite temperature. In the phase diagram of this model there is the so-called 
Nishimori line which intersects the phase boundary at a multicritical point. It is known that the 
correlation functions obey many exact identities on this line. We use a supersymmetry method to 
treat the disorder. In this approach the transfer matrices of the model on the Nishimori line have an 
enhanced supersymmetry osp(2n -I- 1 1 2n), in contrast to the rest of the phase diagram, where the 
symmetry is osp(2n | 2n) (where n is an arbitrary positive integer). An anisotropic limit of the model 
leads to a one-dimensional quantum Hamiltonian describing a chain of interacting superspins, which 
are irreducible representations of the osp(2n + 1 | 2n) superalgebra. By generalizing this superspin 
chain, we embed it into a wider class of models. These include other models that have been studied 
previously in one and two dimensions. We suggest that the multicritical behavior in two dimensions 
of a class of these generalized models (possibly not including the multicritical point in the RBIM 
itself) may be governed by a single fixed point, at which the supersymmetry is enhanced still further 
to osp(2n-|-2 I 2n). This suggestion is supported by a calculation of the renormalization-group flows 
for the corresponding nonlinear sigma models at weak coupling. 



I. INTRODUCTION 

For many decades Ising models served as the simplest 
nontrivial models for the description of magnetically or- 
dered phases and phase transitions between them. This is 
true both for pure models and for Ising models with ran- 
domness. In particular, in the context of the spin glass 
problemH the relevant Ising models have random bonds 
of both signs (ferro and antiferromagnetic) . This leads 
to frustration and the possibility of spin glass order. 

In this paper we consider a classical random-bond Ising 
model (RBIM) of Ising spins Si — ±1 on the two- 
dimensional (2D) square lattice with the Hamiltonian 
(/? = 1/T is the inverse temperature) 



KijSiSj, 



(1.1) 



where the bold indices i = {ix,iy) and j = {jx,jy) de- 
note 2D vectors of integer coordinates of the sites of the 
lattice, the summation is over distinct nearest-neighbor 
bonds (i.e. pairs), and the coupling constants Kij are in- 
dependent random variables drawn from the distribution 

P[Kij] = (1 -p)<5(ifij - K)+pd{Kij+K). (1.2) 

In words, the couplings Kij are ferromagnetic {K > 0) 
with probability 1—p and antiferromagnetic with proba- 
bility p. Notice that K varies inversely with T. In what 
follows we will occasionally also consider Ising models 
with other distributions of the bond strengths. For sim- 
plicity, in most cases, where it cannot lead to confusion, 
we w ill s imply call the model with the binary distribu- 
tion (1.2) "the RBIM". Later, we will also consider the 



anisotropic generalization of the model, in which K takes 
different values on bonds in the x and y directions. 

Let us summarize some of what is known about this 
model. The phase diagram of this model is still somewh; 
controversial, but is widely believed to be as in Fig. npTl 
First we note that for p = 1, we have a pure antifer- 
romagnetic Ising model, which can be mapped onto the 
ferromagnetic case by sending Si — > —Si for i on one 
sublattice. More generally, this transformation is equiv- 
alent to sending p ^ ^ ^ P- Hence we need show only 
the region < p < 1/2. The solid line is a phase bound- 
ary which separates the ferromagnetically-ordered from 
the paramagnetic phase. Fig. |l| can also be viewed as a 
schematic renormalization group (RG) flow diagram, in 
which the intersection points labeled Tc (corresponding 
to Kc = 0.44 . . ., the pure Ising transition), pc — 0.12 (the 
T = transition) , and another point N are viewed as RG 
fixed points that govern the critical behavior for the por- 
tions of the phase boundary shown as flowing into these 
points (N is an unstable, hence multicritical, point). In 
2D, it is generally believed that no spin glass phase ex- 
ists at finite temperature. At zero temperature, long- 
range spin-glass (Edwards- Anderson) order exists triv- 
ially when the distribution of bonds is continuous, since 
in a finite system there is, with probability one, a unique 
ground state, up to a reversal of all the spins. (Taking the 
thermodynamic limit in a fixed sample is a very subtle 
problem; for a very recent discussion, see Ref. ^ and ref- 
erences therein.) However, for the discrete distribution 
with bonds taking values K, —K, assumed here, the ex- 
istence of such order in the region p > pc '\s not clear, 
because there will be many degenerate ground states. 
There is evidence for power-law spin-glass correlations 
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FI G. 1 . Phase diagram of the random-bond Ising model of 
Eq. (1.1), in terms of T oc 1/K and p. 



at r = in this regioiMtll. In three or more dimen- 
sions, there is a spin-glass ordered phase at temperatures 
below some temperature TsgIp); and which extends up 
to p = 1/2. The three phase boundaries meet in a mul- 
ticritical point at some T < Tc, p < 1/2. The point N in 
2D is in some sense a remnant of this multicriticality in 
higher dimensions. 

An interesting feature of this phase diagram is the so- 
called Nishimori line (NL), shown dashed. Such a line 
can be defined for a broad class of distributions PlKa], 



and in our special case (1.2) is given by the condition 



l-2p = tanhi^. (1.3) 

Nishimorfl found that the model on the NL has a local 
( "gauge" ) Z2 symmetry. Using the symmetry he showed 
that the internal energy of the RBIM is analytic along the 
NL, and also established a special case of the following 
identities for correlation functions (proven generally in 
Ref. 0) which hold on the NL 
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(1.4) 



Here {Si-, , . • . , Si, } is any set of the Ising spins, the angu- 
lar brackets denote the thermal average, and the square 
brackets the average over the distribution of bonds -P[/^]. 
Georges et aM showed that Nishimori's result concerning 
the internal energy may be rederived using a supersym- 
metric formulation. Nishimori later argued, and Kitatani 
showedQ, that the ferro-para phase boundary is vertical 
below the NL. This is supported by the later numerical 
works among Refs. ||,^. „ 

Le Doussal, Georges and Harriall established that in 
any dimension the NL goes through the multicritical 



point A'' on the ferro-para phase boundary and is one of 
the exact RG trajectories near this point, the other tra- 
jectory being the vertical tangent to the phase boundary 
at N. The RG eigenvalues along these trajectories and 
the corresponding, critical exponents were estimated by 
Singh and AdlerE3 from high-temperature expansions of 
high order. Two of their exponents in 2D are very close 
to those of classical 2D percolation, as are the results 
in Ref. ||. We emphasize that, while Nishimori's results 
do not apply to every possible distribution of disorder in 
the RBIM, the implications for universal critical proper- 
ties must hold throughout a universality class, and it is 
believed that this class includes the generic RBIM mul- 
ticritical pointQ. 

More recently, Cho and FisherEj proposed a network 
model similar to the Chalker-Coddington model used 
to desciibe transitions between integer quantum Hall 
plateaual^. The Cho-Fisher model was supposed to be 
in the same universality class as the 2D RBIM. They 
simulated their network numerically and found a phase 
diagram with a multicritical point, with numerical crit- 
ical Gjoponents close r-to the ones found by Singh and 
AdleiO. In her thesisE-3, Cho also simulated another net- 
work model which is precisely equivalent to the RBIM, 
and again found exponents close to some of those in Ref. 

0- . . 

Here we will briefly sketch a version of the argument 
of Ref. 1^, since it provided motivation for our work. Let 
us consider the random-bond Ising model with Gaussian 
disorder, 



P[Kij] - 
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V2^ 



exp[-{Kij^Kof/{2A% (1.5) 



so the mean of Kij is Kq, and the standard deviation is 
A. Taking the partition function 



Z = ^ expy^KijSiSj, 

{s.} (ij) 

we replicate and average to obtain 



(1.6) 



Ko ^ StS^ + Y1 StS^StS'^ 



{ij),ab 



(1.7) 



where a, 5 = 1, . . . , n. For finite n this has the form 
of the Ashkin- Teller model, consisting of n coupled Ising 
models. Now compare this with the replicated spin-glass 
model, which is obtained by setting Kq = 0: 



[Z™] = ^ exp 



^A^ J2 StS^SlS 



{ij),ab 



(1.8) 



with a, & = 0, . . . , m— 1. This model has a gauge symme- 
try: it is invariant under site-dependent transformations 
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5° —S^ for all a and any set of Vs. This local Z2 gauge 
symmetry can be fixed by setting all 5'? = 1 for all i, for 
one value of a, say a = 0. Then if m = n + 1, we obtain 
the random-bond partition function with Kq = , up to 
constants. On this line, which is the NL for the Gaussian 
case, there are still consequences of the underlying gauge 
symmetry. The m — n + 1 replicas are still on an equal 
footing, and this leads to the identities ( |l.4| ) on the NL, 
as follows. In the gauge-unfixed model, only correlation 
functions containing an even number of S'^'s at each site 
are nonzero. The correlation functions are invariant un- 
der permutations of the replicas, and so independent of 
whether or not a = is among the components. Since 
different replica components represent distinct thermal 
aver ages in the n — > limit, we obtain the identities 
(1.4), and others, on gauge fixing. The Z2 local gauge 
symmetry in the replica formalism should not be con- 
fused with that of Nishimori, who did not use replicas; it 
is the enlarged permutational symmetry of the replicas 
that corresponds to Nishimori's arguments. Off the NL, 
the identities are lost, but the model can still be written 
as a gauge-fixed version of a system with n + 1 replicas 
and a local Z2 gauge symmetry. 

The preceding argument shows that, in the replica for- 
malism, the NL is special because it possesses a larger 
permutation symmetry S„+i in place of the usual S„. As 
we saw, even off the NL, an additional "zeroth" replica 
spin can be introduced into the model, along with a gauge 
symmetry that can be used to remove the unwanted de- 
grees of freedom, and further on the NL line the zeroth 
spin is symmetric with the others. Now from work ex- 
tending back to Onsager, the 2D Ising model can be 
written in terms of free fermions, which become Majo- 
rana (real Dirac) fermions in the continuum limit, and 
furthermore Ashkin- Teller models can be represented rby 
interacting Majorana fermions, with 0{ti) symmetryll3. 
Hence we are led to conjecture that, in the replicated 
fermion representation, it is possible to introduce an ad- 
ditional "zeroth" fermion, together with a local Z2 gauge 
symmetry to remove the unwanted degrees of freedom, 
and that on the NL, we should find a larger 0{n + 1) 
symmetry. In this paper we demonstrate that this in- 
deed occurs, though we take a different route to do so. 
We consider the binary distribution above, and we use 
supersymmetry rather than replicas, so no limit n — > 
need be taken. However, the corresponding result for 
replicas is contained in our results. The network models 
such as the Cho-Fisher model also depend on the fermion 
representation of the Ising model, and so we can also con- 
sider these models in our framework. We find that the 
models can be viewed as supersymmetric vertex models, 
or by using the anisotropic limit as Hamiltonian chains, 
which act in irreducible representations of the relevant 
symmetry (super-)group (which is enlarged on the NL), 
and thus are quantum spin chains, and possibly can also 
be viewed as the strong-coupling region of a nonlinear 
sigma model. This greatly enhances the similarity of the 
problem to the integer quantum Hall effect transition, 



and to other random fermion problems. However, we 
find that the NL does not fall into a recent list of non- 
linear sigma models that cacrespond to random matrix 
ensembles in such problemsED. Our results also apply to 
certain one-dimensional fermion problems. 

In random fermion problems, including those arising in 
disordered superconductors, it is usual to attempt clas- 
sifications, based on symmetries, for generic probability 
distributions, as in Re f. O for random matrices. The en- 
sembles found in Ref. |17| for disordered superconductors 
differ from the standard ensembles because of the lack of 
a conserved particle number, and because zero energy is a 
special point in the spectrum. The second-quantized non- 
interacting quasiparticle description can in each case be 
replaced by a "first-quantized" formalism involving a sin- 
gle matrix, which must satisfy certain discrete symmetry 
and symmetry-like conditions, which distinguish the en- 
sembles. One such class, termed class D in Ref. is for 
problems with broken time reversal and no spin-rotation 
symmetry. This corresponds to the symmetries of the 
fermion representation of the RBIM. In this class, a non- 
linear sigma model analysis in two dimensions indicates 
that there is a metallicj^hase in which the fecmion eigen- 
states are extendedliaO. Senthil and Fisherlij discussed 
a scenario in which such a phase occurs in the RBIM, 
as an intermediate phase between the para- and ferro- 
magnetic phases, at low T in the region labeled "Para" 
in Fig. |l|, and with its bordering phase boundaries (one 
of which is the low T phase boundary shown) meeting 
at the multicritical point N on the NL. (This is the re- 
gion sometimes claimed to be some kind of spin-glass- 
like phase in the RBIM literature.) They suggested that 
the phase would be characterized by the absence of long 
range order in the mean of either the ferrcapagnetic Ising 
or the dual disorder variable correlationaij. Presumably 
such decay would also hold for the mean square (spin- 
glass) correlations. It is not clear if this is consistent 
with the T — > analysis of Ref. |ll|. An alternative sce- 
nario is that a finer analysis of problems with broken 
time-reversal symmetry is needed, and that the nonlin- 
ear sigma model appropriate for class-D does not apply 
to the RBIM. Indeed, a recent paperEil emphasizes that 
the target manifold of the class D model is not connected, 
and that consequently there can be domains of the two 
components or "phases". It follows that additional pa- 
rameters are required in order to fully parametrize the 
systems, independent of those familiar for sigma models 
with connected targets. This allows for a much richer 
phase diagram and transitions in this symmetry class. 
We show in this paper that this is connected with the 
structure we uncover in the random-bond Ising and net- 
work models. In another paperE3, it is argued that the 
metallic phase cannot occur in a RBIM with real Ising 
couplings. 

Problems of random noninteracting fermions are 
among the better understood of disordered systems, and 
while many results are numerical, in some cases there are 
even exact results for critical properties in 2D. By mak- 
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ing contact between the RBIM and other random fermion 
problems, and casting them in a common language, we 
hope to gain understanding of this disordered classical 
spin problem. At the same time, the RBIM provides 
an example that may shed light on previously-unknown 
classes or ensembles of random fermions. The analysis 
presented in this paper does not resolve all aspects of the 
broken-time-reversal symmetry class, but it does show 
that the NL is a special subclass. 

We now give an outline of the paper. Sections ||, 
III, IV and VI contain the main technical work. They 
show how a fermion representationc2l can be used for 
the RBIM, and how bosons are also introduced to can- 
cel the iii\ferse partition functions, via supersymmetry 
(SUSY)E3~e1I. The bosons live in a space with an indef- 
inite metric, a common feature of SUSY methods (see 



for any sign of Kij . For positive Kij = K > we denote 



Refs. E6 31 



and compare Refs. |28|-|30D. On the NL, a 
larger SUSY algebra is found. As an application of this 
enhanced SUSY, we use it in Appendix ^ to rederive the 
equalities ( |l.4| ). An anisotropic limit of the model relates 
it to a Hamiltonian for "superspins" on "split" sites, two 
for each original site. The Hamiltonian on the NL can be 
generalized in a natural way, and we introduce a phase 
diagram for the generalized models; on one line in this 
diagram the Hamiltonian posspses pS- still larger SUSY. 
Section discusses a ID modelEaHn^, and shows that it 
poss esse s a single "Nishimori" point of higher SUSY. In 
Sec. |VIl| we consider the model of Cho and FisherEj, and 
show that it does not correspond precisely to the RBIM, 
and nor does it have the larger SUSY of the NL. In Sec. 
VIII, we argue that much, or possibly all, of the critical 
surface in the phase diagram in the higher SUSY gener- 
alized models is in the universality class of the point with 
even higher SUSY mentioned above. This is supported 
by consideration of nonlinear sigma models that corre- 
spond to the spin chains, and a weak-coupling calculation 
for these shows a renormalization-group flow towards the 
higher SUSY theory. We comment on general network 
models and nonlinear sigma models with the symmetries 
of class D. Appendix ^ gives details of a representation 
of SUSY that we use. Appendix |c| contains the details 
of the calculation of the beta functions at weak coupling 
in the nonlinear sigma model. 



II. TRANSFER MATRICES AND 
SUPERSYMMETRY 

In this Section we will express the Ising model transfer 
matrices in terms of fermionic replicas, and then intro- 
duce bosons to make the system supersymmetric. This 
allows us to consider (in Sec. [II) averages over quenched 
disorder without taking the replica n — > limit. 

First, we set up some notation. We define the dual 
coupling K by 



tanh Ki 



(2.1) 



Kij — K* , and for negative Ki 



-K < 0, Kij = K* 



i7r/2. We assume free (not periodic) boundary conditions 
on the Ising spins in the horizontal [x) direction, and 
periodic in the vertical {y) direction. As is well known 
(see, for example, Ref. |36|), the partition function of the 
nearest-neighbor Ising model may be written as the trace 
of a product of row transfer matrices: 



Z = Tr W_T^{iy)Th{iy 



(2.2) 



Here iy is an integer coordinate of a row of sites. The 
row transfer matrices, T^{iy) for vertical, and Th{iy) for 
horizontal bonds, d o no t commute with each other, so 
the product in Eq. {2fl) must be ordered such that the 



row coordinate iy increases from right to left. The row 
transfer matrices may in turn be written as products of 
the transfer matrices for single bonds: 



(2.3) 



The T^i's for different and the same iy commute, and 
similarly for the T/ii's. The trace represents the periodic 
boundary condition in the y direction. 

Following Ref. ^ we write the vertical transfer matrix 
for a single vertical bond between the Ising spins at points 
i and i -)- y as 



T 

J. 7). 



cosh K, 



i,i+y 



- exp {ki^i+yaf^ , 



(2.4) 



where the "light" index i denotes a site on the ID lattice 
corresponding to the vertical row containing the original 
site i, and cr^, cr^ and are Pauli matrices. Similarly, 
the horizontal transfer matrix for a single horizontal bond 
between the Ising spins at points i and i -I- a; is 



Thi = exp (i^i,i+£CTf CT^+i) 



(2.5) 



Note that before the averaging over the randomness the 
transfer matrices explicitly depend on the corresponding 
bond, and therefore are labeled by the bold 2D indices. 

The transfer matrices act in tensor products of two- 
dimensional spaces at each horizontal coordinate ix- 
These 2D spaces may be realized as Fock spaces of 
fermions on a ID chain of sites. This fermionization is 
implemented by the Jordan- Wigner transformation re- 
lating Pauli matrices to fermionic operators. To use this 
transformation we first make a canonical transformation 



The Jordan- Wigner transformation reads 

- 24c^ - 1, 



(2.6) 

(2.7) 
(2.8) 
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where c] and Ci are canonical creation and annihilation 
fermionic operators. In terms of these operators the 
transfer matrices for individual bonds become 



T,, 



cosh K, 



cxp -2Xi,i+y(c|q-l/2) , (2.9) 



= exp (ifi,i+i(c| ~ cj(c|+i + q+i)) . (2.10) 

Next we replicate the fcrmions. The number of replicas 
has to be even, because of the symmetry of the bosons 
to be introduced below, so we denote it as 2n. We label 
the replicas by Greek letters. The replicated transfer 
matrices become 



T„i = (2 cosh if)2" exp (^-2Ki^i+yNF,'^ , (2.11) 
Thi = exp (2Ki^i+^XF^) , (2.12) 



(2.13) 



where we defined 



2n 

Npi = rial, 

a=l 
2n 

a=l 



1-ai — C^^C^i, 



(2.14) 



(2.15) 



The quadratic forms Npi and Xpi are invariant un- 
der the orthogonal transformations mixing the fermions, 
which becomes especially transparent if we introduce two 
sets of 2n real fermions per site as 



V2 



(2.16) 



V2i ' 
These fermions satisfy 

{Vat, Vflj} = {^ai,^f3j} = S.jSap, {Vai, C/3j } = 0. (2.17) 

Terminologically, we note here that any set of self-adjoint 
operators, say ipa, a = 1, . . . M, for some M, with anti- 
commutation relations {ipa-, ^b} = Sab, constitutes a Clif- 
ford algebra. For us, the set of ^'s, either for one or for 
many sites, or similarly of ?7's, or a combination of these, 
are all Clifford algebras. A little of the general theory 
of these algebras will be used later. In terms of these 
fermions, or Clifford algebra generators, the quadratic 
forms become 

Npi = iriaiiai + n, Xpi = iriaiCa,i+i, (2.18) 

where from now on we assume that repeated indices from 
the beginning of the Greek alphabet (a, f3, etc.) are 
summed from 1 to 2n, unless stated otherwise. 

The generators of the global symmetry algebra so(2n), 
in this notation, are Y^iiVaiVpi + iai£,i3i), for pairs a, /?, 
and because of the anticommutation relations we may 



take only a < P, corresponding to the antisymmet- 
ric 2n X 2n matrices. These generators commute with 
Npi and Xpi, proving that the transfer matrices are in- 
variant under so(2?i). The replicated partition function 
Z^", which is now given by a trace in the 2ri-component 
fermion Fock space, is invariant under so(2n). Note that 
we capitalize the name of the group or supergroup, such 
as S0(2n), but not the name of the corresponding Lie 
(super-)algebra, such as so(2n). 

The supersymmetric counterpart of the fermionic al- 
gebra so(2n) is the symplectic algebra sp(2n). This mo- 
tivates the introduction of bosons with this symplectic 
symmetry as follows. We start with 2n complex "sym- 
plectic" bosonic operators Sai (and their adjoints sl^^), 
satisfying 



(2.19) 



where Jap is a non-singular real antisymmetric 2n x 2n 
matrix. (It is because the number of bosons must be even 
that the number of fermions must be also.) Without loss 
of generality (by appropriate change of basis) this matrix 
may be taken to be (in block form) 



J 



1„ 
-1„ 



(2.20) 



where 1„ is the n x n identity matrix. We also need 
to define the vacuum state for our bosons. We will see 
that the SUSY requirement makes this choice essentially 
unique, and the resulting space of states has indefinite 
metric (some states have negative squared norms). 

The bosonic counterparts of the forms Npi and Xpi 
are the symplectic forms 



(2.21) 
(2.22) 



To parallel the fermionic case we also introduce two sets 
of 2n real bosons per site as 



9q 



V2 

These bosons satisfy 



V2i 



(2.23) 



ikjJap, [qai,ri3.j\^Q. (2.24) 



(These have the form of the commutation relations for 
canonically conjugate coordinates and momenta.) In 
terms of the real bosons the forms (2.22) are 



Bi 



X 



Bi 



' Oil Jai3qp,i+i- (2.25) 



The generators of the global symplectic symmetry al- 
gebra sp(2n) are J^iilaiQ/Si + fai'<'pi), where because of 
the commutation relations we may use only a < (3, cor- 
responding to symmetric matrices. These operators are 
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the generators of global linear canonical transformations 
on the g's and r's. The forms Nsi and Xsi, and hence 
the transfer matrices, are invariant under this algebra. 

We now address the question of the bosonic vacuum. 
We will find it using the requirement that the spectrum 
of the bosonic form Nsi is the SUSY counterpart of the 
integer spectrum of the fermionic form Npi- In this case 
we will have 



STr exp 



1. 



(2.26) 



This condition is essential in the SUSY approach to en- 
sure that the partition function of the RBIM is unity for 
any realization of the disorder. Here we used the no- 
tation STr for the supertrace in the space of states of 
our problem. We will now discuss how this supertrace is 
defined. 

In general, the supertrace in a super-vector space must 
be defined using the notion of a grading for the states 
(or vectors). This can be done by choosing a basis and 
then defining one subset of basis vectors as "even" , and 
the remainder as "odd" , vectors. The vector space then 
contains two complementary subspaces of even and odd 
vectors, respectively; the zero vector, and linear combi- 
nations of vectors from both subspaces, are viewed as 
having no definite grading. The vectors are then said to 
be Z2-graded. Operators on the vector space can like- 
wise be classified as even or odd, according to whether 
they preserve or reverse the grading of basis vectors on 
which they act; usually, only operators for which this rule 
gives a consistent answer (those with well-defined grad- 
ing) are of interest. Thus the grading is usually treated 
as a superselection rule. The supertrace STr Y of an even 
operator Y is then defined, like an ordinary trace, as the 
sum of the diagonal matrix elements of K in a basis of 
even and odd vectors, except that for the supertrace, the 
diagonal elements in the odd basis vectors are weighted 
by a minus sign. (Note that the matrix elements Yij of 
an operator Y are obtained as the coefficients in the sys- 
tem of equations Y\J) = '^jYij\I), where |/), / = 1, 2, 
. . . , are the basis vectors, without using an inner prod- 
uct on the vector space.) The supertrace has a number 
of nice properties, like the ordinary trace in an ordinary 
vector space; in particular a form of the cyclic property 
stih holds, STr AB = ±STtBA, with a -h if both A 
and B are even, and a — if both are odd operators. The 
grading and the supertrace are needed in connection with 
supersymmetry algebras, but otherwise do not necessar- 
ily have to be considered. We also note here that it is 
possible to form a graded tensor product of graded vector 
spaces, in a way that preserves the grading. 

In this paper, most of our constructions use a Fock 
space. In a Fock space generated by boson and fermion 
operators acting on a vacuum, there is a natural grading, 
defined using an occupation number basis, in which states 
are even or odd according as the total number of fermions 
(of all types) is even or odd. However, we will not use 
such a grading to define the supertrace above. The reason 



is that we have already introduced the ordinary trace in 
writing the partition function for fermionic replicas; in 
this trace, all diagonal matrix elements are taken with 
weight -|-1, including those in states with an odd num- 
ber of fermions. It is of course quite standard to use an 
ordinary trace even when dealing with fermions, which 
have a natural grading. The natural grading is used in 
defining a tensor product, such that fermion operators 
on different sites (i.e. in different factors in the tensor 
product) anticommute. These are the tensor products 
usually used by physicists for second quantized fermion 
problems. Each time we write a tensor product of spaces, 
it will be the graded tensor product using the natural 
grading that we mean. There is nothing wrong with the 
use of the trace, unless we are concerned about SUSY. 
The grading that we use in introducing SUSY into our 
representation is defined by specifying that states with 
an even (odd) number of bosons are even (odd), and so 



STr 



= Tr(-l)S' 



(2.27) 



For states with no bosons, this reduces to the usual trace. 

Now that we have defined the supertrace, we must ar- 
range to satisfy the condition (2.26). The form Nbi may 
be diagonalized with the transformation to two other sets 
of n complex bosons. Namely, we define 



^Sp^^rtji 



V2 ' V2 



(2.28) 



and the adjoint operators, where the index fj. (and other 
indices from the middle of the Greek alphabet, like i^, 
etc.) runs from 1 to n. These bosons satisfy 

Ki' ^^i^j] = '^y'^^'^' t^M*' " (2-29) 

and the rest of commutators vanish. In terms of these 
bosons we have 



(2.30) 



If we introduce the vacuum for a and a bosons in the 
usual manner 



(2.31) 



then the spectrum of Nbi is the non-negative integers, 
which is the SUSY co unterp art of the spectrum of Npi. 
This ensures that Eq. (|2.26| ) holds. But the price to pay 
is that the states with an odd number of a bosons have 
negative norms. By another choice of the vacuum we 
could avoid negative norms, but then we would not have 
the supersymmetry. With these definitions, we have now 
defined a Fock space J-', which is a tensor product of 
Fock spaces at each site, J- = ®i!Fi in an obvious nota- 
tion. The tensor product of Fock spaces is defined using 
the natural grading, however our choice of grading also 
behaves well in the product; the grading of states is de- 
termined by the product of the "degrees" (= ±1 for even. 
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odd respectively) of the states on the sites, because boson 
numbers add. Note that fermion operators are viewed as 
even in our grading. 

The transfer matrices including fcrmions and bosons 
supersymmetrically are now 



r„i = cxp y-2Ki^i+yNstj , 

Thi - exp (2Ki^i+iXs^) , (2.32) 

where the subscript S stands for "supersymmetric" , 

Ns^ = Nf^ + Nb^. Xs^ = Xf^ + Xb^. (2.33) 

The SUSY transfer matrices are invariant under the or- 
thosymplectic superalgebra osp(2n|2n), since the forms 
Nsi and Xsi commute with the generators J^ii^ai^Pi + 

VatVPt), J2iilai(l0i + 'TaiTpi), and X]j(Cai'?/3i + Vatrpi) of 

osp(2n|2rt). The last set of generators are the "odd" 
(with respect to either grading), fermionic, or supergen- 
erators of the superalgebra, and a and (3 can take ar- 
bitrary values there. Note that, in a superalgebra, two 
even operators obey commutation relations, two odd op- 
erators anticommutation relations, and an even with an 
odd generator obeys a commutation relation. Thus the 
definition of the superalgebra structure again involves the 
grading. The definition of the supertrace also respects 

supersymmetry. 

The condition ( 2.26| ) applies when only vertical-bond 
transfer matrices are present. To prove the supersym- 
metry of the full problem, namely, that the supersym- 
metrized partition function ZgusY (the supertrace of the 
product of supcrsymmctrized transfer matrices) is unity 
for any realization of the disorder, we use a graphical rep- 
resentation of Z . Imagine a high-temperature expansion 
of the horizontal transfer matrix Tm, where we expand it 
in powers of i^i.i+x • At each horizontal row the operator 

2Xsi = (c^ - c„J(cJ, + c„_,+i) 

may create or destroy two particles at neighboring sites, 
or it may transfer one particle between neighboring sites. 
The vertical transfer matrices are diagonal, so they only 
propagate particles in the vertical direction. 

We can represent these processes graphically by lines 
starting and ending at lattice sites (creation and anni- 
hilation), joining the horizontal pairs of sites (hopping), 
and joining vertical pairs of site (propagation due to T^). 
When two such lines start at a pair of neighboring sites 
on a row and end at another pair of neighboring sites on 
a row, we call this a closed loop (similar to closed loops 
in high temperature expansion of the pure Ising model). 
Then the supersymmetrized partition function is equal to 
1 plus the sum of contributions of all closed loops. Then 
we have to prove that for a closed loop the contributions 
of fermions and bosons cancel each other. 



Let us take the smallest possible loop, where two parti- 
cles are created and destroyed on two adjacent rows. This 
is represented by two short vertical lines between two 
neighboring pairs of lattice sites. For a given fermionic 
replica, say 1, this loop contributes the following term: 



The corresponding bosonic contribution is 



(2.35) 



(2.36) 



To evaluate this expression we note the following. From 
the definition of the symplectic bosons it follows that 



-<5.,4.. (2.37) 



We use these relations to pull the expo nentials through 
to the vacuum on the right in Eq. (2.36), which becomes 
then 



i^n-|-l,i+ 



i|0)e 



+ it,i + ii! + v) 



(2.38) 



Next we notice that, as a consequence of Eq. (2.28) and 
the definition (2.31), the vacuum state |0) is annihilated 
by s^^. Then we commute the operators in the first fac- 
tor in the last expression, after which it becomes exactly 
opposite to the fermionic contribution (2.35). This argu- 
ment is easily generalized to arbitrary loops (including 
those that wrap around the system, thanks to the defi- 
nition of the supertrace), and proves that the supersym- 
metrized partition function is indeed equal to one, for 
any realization of the disorder. 



III. AVERAGING AND ENHANCED 
SUPERSYMMETRY ON THE NISHIMORI LINE 

In this Section we perform the average with respect to 
the distribution P[K] and find that on the NL the av- 
eraged transfer matrices have enhanced supersymmetry 
osp(2n + 1 1 2n). In Appendix [a| we use this enhanced 
SUSY to rederive the equality (1.4) for the Ising correla- 
tors. 

Here we need some more notation. We introduce a 
parameter of the form of the Ising coupling, L, and its 
dual L* , related to the probability p: 



1 — 2p = tanh L 



-2L' 



(3.1) 



In terms of L the equation ( |l.3| ) of the NL \s L = K. 
Below the NL i < iiT, and above the L> K. 



7 



Since the couplings Kij are independent, we can aver- 
age transfer matrices for different bonds separately. For 
a vertical transfer matrix this gives (recall that the dis- 
order averages are denoted by square brackets) 

Tu = [Tyi] - exp (-2if*7Vs,) (l - p + pi-lf'' 

(3.2) 

Note that after the averaging the translational invariance 
is restored, and this allows us to label the average transfer 
matrices by ID ("light") indices. 



The value of the last factor in Eq. (3^) depends on the 
value of Nsi in the state, on which Tu acts. For an even 
Nsi it equals 1, for an odd Nsi it give s 1 — 2p = e~^^ . 
Then we can rewrite the operator ( |3.2D in a slightly dif- 
ferent form. Namely, we introduce additional (zeroth) 



fermionic state and operators 



and consider the 



subspace J-^ given by the following constraint: 

^Si = ^Si + noi = even, 
where noi is defined in analogy with Uai : 



(3.3) 



That is, the number of fermions plus bosons on each site 
must be even. There is a one-to-one correspondence be- 
tween the states in this subspace and the original Fock 
space !Fi. This correspondence is illustrated for the case 
n = 1 in the Table |. The grading in the space T' — ®iJ-[ 
is taken to be the same as that in JF, which was not the 
natural grading. However, we see tha t in a Fock space 
with a constraint of the form of Eq. (|3.3| ), the number 
of bosons is odd if and only if the number of fermions 
is odd (this is true for each site and also for the tensor 
product). Hence, our grading on T' is the same as the 
one obtained from the natural grading on the larger Fock 
space with 2n+l fermion species, when restricted to the 
subspace. Now we can replace Tu by the operator 



Ti^ = exp (-2K*Nsr - 2L*no, 



(3.5) 



which has the same matrix elements in the constrained 
subspace JF^', as Tu had in the original space (be- 
tween the corresponding states). From now on in this 
Section we will denote transfer operators acting in the 
constrained spaces J^' by a prime. 

For a horizontal transfer matrix the averaging gives 



Toi — \Thi\ — 



cosh{2KXs^ + L) 
coshL 



(3.6) 



To find the corresponding operator in JF' we need to es- 
tablish some substitution rules for basic operators. 

Single creation and annihilation operators like c 
which are quite legitimate in the space .Fi, do not act 
within T[. Using the correspondence between the states. 



t 



|0> 

4lo) 
cL4lo> 
«L|o> 



■-2 



|0) 

4lo) 
4lo) 
«L|o) 
«Llo) 
«L|o) 

4aL|0) 



TABLE I. Correspondence between states in the spaces JF^ 
and T'. 



given in Table |, it is easy to establish, that in JF/ the op- 
erator cjj must be replaced by the operator (cj^ — Cq^)c\^. 
However, this operator is bosonic (i.e. even with respect 
to the natural grading on J^'), so in the tensor product 
.F' = ®iJ^l it will not have the anticommutation proper- 
ties that cl, had in To make it fermionic in the total 



(3.4) space !F', we need to attach to it a string: 



(3.7) 
(3.8) 



Repeating this argument for other operators, we obtain 
the following rules of substitution: 



- coi)4»Soi, 

(4 ~ CqJc^i^Oj, 

(4 + Coi)«IiSo., 

(4 + Coi)V^O., 

(4 + coi)4iSo., 

(4 + Coj)«uiSo.- 



(3.9) 



As an alternative to looking at all the states, the cor- 
respondence can be established by verifying that the 
right-hand sides of these expressions have the same (anti- 
)commutators as the left-hand sides. It follows from these 
rules that any product of an even number of creation 
and/or annihilation operators in J^i remains the same, 

when going to J-^. 

Now note what happens with Xai (see Eq. ( 2.15| )) upon 
transition from J- — to T': 



2xc\ 



i)(4/j+l 



l) 



= (4. ^ 

-(4-co.)(4.-c„.)(-i)"°"+^ 

^('^0,1-1-1 ~ '^0,1-1-1 ) (4,1-1-1 + '^a,i+l) 
= AxQiXai^ (3.10) 
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where xoi is defined similarly to Eq. (2.15) 



^Oi = - CoJ(cJ ^+i + Co^,+i) = im^^O,^+l■ (3.11) 

With these substitution rules established, we can see 
that in the space J-'' the operator corresponding to T2i is 
given by 



T' 



Oi) 



coshi 



(3.12) 



which easily follows from the substitution rule_( j.lfl) an d 
the fact that a;^. = a;^. = i/4_ pj-om Eqs. (^), {^.12\j, 



we see that t he tr ansfer matrices commute with the con- 
straint, Eq. (3.3). In the full space of states that in- 
cludes states of odd, as well as even, fermion plus boson 
number at each site, there are local Z2 operations given 
by (— l)^si. This is therefore a gauge symmetry under 
which the allowed states and transfer operators (and also 
physical obse rvab le s) mu st be invariant. 

The forms (3.5), ( 3.12| ) are very convenient for the dis- 
cussion of the symmetry properties of our model on the 
NL. Indeed, we see that on the NL, where K = L, the 
operators T/ become 



T' 



T' 



exp (-2K*N's,^, 



cos\y{2KX'^ 



Sil 



coshiiT 



(3.13) 
(3.14) 



where 



X'g^ — Xsi + Xoi = ir]ai(.a,i+l — raiJaf3qfj,i+l, (3.15) 

and the Latin subscripts from the beginning of the al- 
phabet denote the fermionic indices running form to 
2n. 

On the NL, the expressions ( 3.13 ), ( 3.14 ) have en- 
hanced supersymmetry: they are now invariant under 
an osp(2n -I- 1 1 2n) algebra. The generators of this al- 
gebra have a similar form as before, but involve the 
2n-|-l fermion operators: J^zi^at^bi+VaiVbi), Y.i{la^1l3i + 
raiViSi), and J^Maiqpi + VaiTpi)- The last set of gener- 
ators are the odd ones, with respect to our grading, or 
to the natural one on the Fock space of which T' is a 
subspace; we have seen these are equivalent in the con- 
strained subspace. 

As anticipated in the Introduction, this enhanced con- 
tinuous SUSY replaces the gauge symmetry of Nishimori, 
and the enhanced permutational symmetry of the replica 
approach, previousiy known to exist in Ising spin lan- 
guage on the NLqlI. The sy mme try has many conse- 
quences, such as the equalities (1.4) among different cor- 
relation functions on the NL. In appendix we briefly 
show how these equalities may be obtained from the en- 
hanced SUSY exhibited in this section. We have also 
obtained the local Z2 gauge symmetry anticipated in the 
Introduction. 



IV. STRUCTURE OF THE SPACE OF STATES 
AND THE HAMILTONIAN LIMIT 



In this Section we first analyze (in Sec. IV A) the struc- 



ture of the space of states of our quantum problem and 
then take the time continuum limit of our transfer ma- 
trices and obtain a quantum Hamiltonian describing our 
system. This has the form of a spin chain with irreducible 
representations of the sym metry algebra osp(2n -I- 1 1 2n) 
at each site. Then in Sec. IV B we consider a more ex- 



plicit construction of these irreducible representations. 



A. Superspin chain and Hamiltonian limit 

Let us consider the structure of the constrained space 
JFj' ( p.3[ ) , with its natural grading, under transformations 
of osp(2n -I- 1 1 2n) . It is easy to see that JF^' is not ir- 
reducible under this algebra. Rather, it has the struc- 
ture of the tensor product of two irreducible spinors of 
osp(2n-|- 1 1 2n). 

Indeed, let us consider first the fermionic replicas only, 
i.e. the replica approach where n ^ in the end. Then 
we have the modified transfer matrices which are invari- 
ant under the orthogonal algebra so{2n + 1), and the 
subspace they act on is given at each site by the con- 
straint 



Npi + noi = even. 



(4.1) 



This space has dimension 2^", and, under so(2n -f- 1), 
it transforms as the tensor product of two spinors of 
so(2n -I- 1), each of dimension 2". These two spinors can 
be identified as the spaces on which the two parts ^ai^bi, 
VaiVbi of the generators ^ai^bi + TlaiVbi act. The tensor 
product decomposes into irreducible representations of 
so(2n -I- 1) corresponding to each even valu e of N'p- in 
the range to 2n allowed by the constraint (4.1). Simi- 
larly, the orthogonal subspace Np^ = odd is also a tensor 
product of spinors and has a similar decomposition. 

When the bosons are included as in the SUSY ap- 
proach, the two parts qailpi, r^irpi of the sp(2n) gen- 
erators at a site i generate infinite-dimensional spinor 
representations of sp(2n) (sometimes known as meta- 
plectic representations). When the fermions and bosons 
are combined together with the constraint that N'g^ be 
even, the resulting space is a tensor product of irre- 
ducible spinors of osp(2n -|- 1 1 2n). The fact that a single 
such tensor product is involved is the nontrivial part of 
this statement, and is addressed further in Sec. [VB. 
These spinors comprise one lowest-weight representation 
of osp(2n -I- 1 1 2n), which we denote by i?, and one high- 
est weight representation of osp(2n -t- 1 1 2n), which we 
denote by R. Thus we may write T[ = Ri ^ Ri. 

This organization of states suggests a picture of our 
model as a system of "superspins" (spinors R and R of 
osp(2n + 1 1 2n)) sitting in pairs on the sites of the ID 
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h 

-•■ -e- 



h 

—m- -o 



FIG. 2. The graphical representation of the Hamihonian 
( [I.T] ) on the split sites. Superspins in the representations R 
and R are shown as filled and empty circles. The two types 
of coupling, h and k, are indicated. 



lattice. It is convenient to combine the corresponding 
generators of osp(2n + 1 1 2n) into square matrices con- 
sistent with reality properties satisfied by the matrices of 
osp(2n + 1 1 2n) in the defining representation. Namely, 
the generators of osp(2n + 1 1 2n) acting in the represen- 
tation R are combined into the superspin 



G = 



^a£,b — hSab 



(4.2) 



shown here in block ((2n -|- 1) -I- 2n) x ((2n -I- 1) -I- 2n) 
form, and a similar matrix obtained from the generators 
of osp(2n + 1 1 2n) acting in the representation R: 



G = 



VaVb - 



(4.3) 



The generators of the global SUSY of the system are now, 
in matrix form, + Gi). 

Next we take a time-continuum (Hamiltonian) limit of 
the transfer matrices. To do that we actually have to 
start with an anisotropic RBIM, where the vertical and 
horizontal couplings take different values of K. In that 
case we can arrange for the situation when in the vertical 
matrix K* , L* ^ 1, and in the horizontal one K^L 
(we continue to use these notations, so that K* , L* are 
no longer related to K, L by the duality relati on). Then 
we can expand the horizontal transfer matrix ( 3.12 ) as 



T' 



1 + 2K^Xli - 
exp (2K^Xl, 



(4.4) 



Then we can combine all and Tji i^^^o a single "evo- 
lution" operator in imaginary time with the Hamiltonian 



Hs^Y. + A* no.) - fc(X|, + 2\Xs, 



(4.5) 



where we introduced 



h^2K\ k = 2K\ A* = — , A=-. (4.6) 

A A 

The parameters A and A* introduce anisotropy in the 
couplings among the replicas. 



The Hamiltonian (^^) may be rewritten in a very sug- 
gestive form using the superspins G and G: 

Hs=Y. {H^s^ + A*no,) + k str AG,AG,;+i) , (4.7) 

i 

where 



A = diag(A, l4„) 



(4.8) 



is a diagonal matrix representing the anisotropy in the 
superspin space. The supertrace here, denoted str, is 
over the 4n -|- 1-dimensional space as above, with -t- for 
diagonal matrix elements in the 2n+ 1-dimensional block, 
and — for those in the remaining 2n-dimensional block. 
With this definition, the expressions reduce to the ordi- 
nary trace for the replica formalism where only the ^'s 
and ?7's are kept (and then n ^ 0), without an overall 
change in sign. 

Now we should notice that in the anisotropic version of 
the RBIM, there are in general two couplings, K,^, Ky, 
and two parameters for the probabilities, L^, Ly. The 
Nishimori condition becomes two equations, = L^, 
Ky — Ly. Thus the NL is replaced by a two-dimensional 
surface (or 2-surface) in the four-dimensional space, and 
so does not divide the phase diagram into two pieces. We 
will continue to refer to this as the NL. There is presum- 
ably a line on this surface at which a transition occurs. 
The complete phase boundary is three-dimensional, and 
the multicritical behavior is found on a 2-surface on this 
3-surface. The multicritical line on the NL presumably 
lies in the multicritical 2-surface on the phase bound- 
ary. Even though the transfer matrices do not have the 
larger SUSY everywhere on that 2-surface, we presume 
by universality that the higher SUSY fixed point theory, 
to which the multicritical point on the NL flows, controls 
the entire multicritical 2-surface, because anisotropy such 
as we have introduced usually does not affect the univer- 
sality class. 

In any case, on the NL A = A* = 1, and we obtain an 
osp(2n + 1 1 2ri)-invariant (or isotropic) Hamiltonian 



E 



hN'si fcstrG.Gj+i 



(4.9) 



The Eq. (4.9) has the form of a superspin chain with the 
alternating lowest- and highest-weight representations R 
and R, and corresponding superspin operators G and 
G. We can better represent this by splitting the origi- 
nal sites into pairs o f sp lit sites. This is shown in Fig. 
||. The Hamiltonian (O) has two types of couplings on 
the alternating bonds. Both these couplings are antifer- 
romagnetic in nature. This means that the lowest energy 
state for a given bond is the singlet of osp(2n -I- 1 1 2n) 
contained in the decomposition of the tensor product of 
the representations R and R. 
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B. Unconstrained representation of superspins 

The picture of the spht sites carrying irreducible rep- 
resentations is very attractive, but does suffer from one 
difficulty at present. This is that we obtained the rep- 
resentations by introducing an additional zeroth fermion 
coi, cjj, together with a constraint which refers to both 
the split sites that comprise the original site. Here and in 
Appendix ^ we show how to avoid this by use of a differ- 
ent construction. In Sec. VI 
further, introducing a furt 



we will extend this approach 
rer representation which in- 
volves a constraint on each split site. 

The representations R and R can be constructed using 
complex fermionic and bosonic operators without con- 
straints, that is essentially in JF^. For the simplest case 
of the osp(3|2) algebra this is done in Appendix Here 
we note that the complex fermions and bosons used in 
the construction are related to the real ones (apart from 
the zeroth fermion) on the split sites introduced so far in 
the following manner: 



V2 ' 



V2 



V2 



(4.10) 
(4.11) 



All of these operators are canonical, except for the b 
bosons, which are "negative norm" : 



(4.12) 



In terms of these complex bosons and fermion s the 
quadratic forms appearing in the transfer matrices ( 2.32| ) 
look especially uniform: 

^Si = fl^i+lfjii + ffj.iff_i,i+l + + 

Ns^ = fUl + f,J,., + bl^, + b^,b^,. (4.13) 

In this form the subalgebra osp(2n | 2n) on a split site 
is generated by bilinears as before, now of the form 
P, iP?, /V, b^ pb, for R. We saw ear- 
lier that the expressions for Xsi and Nsi are invari- 
ant under this osp(2n | 2n) algebra. However, the ex- 
tension to osp(2ri -I- 1 1 2n) is modified since we do not 
use the zeroth fermion c. Instead, the additional genera- 
tors include string operators such as (—I)"-' , see App. M, 
where expressions for the generators of osp(2n + 1 1 2n) 
in the case n = 1 are given. It is then clear that the 
states on the unsplit sites decompose into a single ten- 
sor product of irreducibles R and R as claimed. The 
string operators again correspond to the difference in 
grading between that natural in and our choice, which 
agrees with the natural one in JF'. With our choice, the 
(anti-)commutation relations obeyed by the generators of 
the larger osp(2ri + 1 1 2n) SUSY are consistent with the 
stated grading, as discussed in more detail in App. 

An advantage of the unconstrained representation of 
the states on the split sites is that it makes the pure 



limit p = transparent. The pure Ising problem in the 
anisotropic time-continuum limit in this representation 
gives a nearest-neighbor "hopping" -type Hamiltonian for 
fermions and bosons, which is a sum over i of Nsi and 
Xsi, with coefficients. This is a lattice version of the 
Dirac fermion and its SUSY partner, the so-called /3-7 
system of bosonic ghosts. However, in the general disor- 
dered case, there are additional terms which we expressed 
previously using the zeroth fermion. In App. |b| we show 
how the GiGi terms in the Hamiltonian can be expressed 
in the present language. The other term, which becomes 
iVcj on the NL, is much more difficult to express in this 



language, and we return to this problem in Sec. VI 



A slight subtlety involved in the definition of the com- 
plex bosons is the following. If we express them in terms 
of the a bosons, which diagonalize the form Nsi, Eq. 
( ^.30 ), we obtain a singular Bogoliubov rotation: 



V2 



-'fit 



V2 



(4.14) 



The singularity of this transformation is seen in the fact 
that the formal expression for the b, b vacuum (on a single 
unsplit site i), defined by b\0) = b\0) = 0, is 



|0) oc exp(— 



|0), 



(4.15) 



where |0) is the a, a vacuum defined in Eq. ( ^.31 ), and 
leads to a series for the squared norm (0|0) which is not 
convergent. A way out of this problem is to regularize 



the Bogoliubov rotation (4.14) as follows: 
bf^i = cos + sin aj^^, 



COSfflO 



(4.16) 



where 4> = i'/4 — w/2 with < w <C 1. With such 
regularized transformation, the b, b vacuum 



1 



cos" 



■exp( 



(-tan0^ |0) (4.17) 



is well-defined and normalized to 1. 

If we now use the regularized relations (4.16), the ex- 
pression for Nsi becomes (to first order in uj) 



Nb, = bl-bl^ + b^^b^i + uj[nu + ni,) - n, 
with bosonic number operators defined as 



''^bi — ~b^fi^^. 



(4.18) 



(4.19) 



The fermionic sector in our formulation is finite di- 
mensional, and there are no similar problems with the 
fermions. However, to maintain the exact cancellations 
between fermions and bosons, we will modify the defi- 
nition of / and / similarly to that for the bosons. One 
effect of this is that the supersymmetric analog of Eq. 
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( [4.17[ ) lacks the factor 1 / co s" cj). Then the osp(2n | 2n)- 
invariant combination ( 4.13| ) takes the form 



(4.20) 



where the fermionic number operators are defined in a 
natural way: 



ifi 



= f'^ f 

•> LilJ Lit ' 



= f-f 



(4.21) 



The term first-order in ui breaks the SUSY down to 
gl(n|n), which is still enough SUSY to ensure cancellation 
of fermions and bosons. 

We can make this term appear more natural by the 
following considerations. It is a regularizer which sup- 
presses contributions to the partition function from high 
fermion and especially boson numbers on any site. We 
can introduce it in a more symmetric way by inserting 
exp uj{nfi + + nt,i + n^J between all the Tn's and 
T2i's in the partition function; to first order in w, the 
effect is the same. Such an insertion is a precaution sim- 
ilar to that often used in network models and nonlinear 
sigma models of localization. The oj term represents a 
non-zero imaginary part of the frequency in those prob- 
lems, and as in the present case breaks the symmetry to 
a subgroup. In the superspin chain language, the opera- 
tor which Lu multiplies is one component of the staggered 
magnetization, the order parameter for the chain. The 
term, with w — > 0, is used just in case this develops a 
spontaneous expectation value, since it picks a direction 
for the ordering in superspin space and cuts off infrared 
divergences. Note that the_ state with each site in the 
vacuum state for the /'s, /'s, 6's, &'s is the Neel state 
corresponding to such order, and is invariant under the 
subalgebra gl(n|n). The symmetry-breaking term will be 
important in Sec. ^ 



V. DIMERIZED LIMIT AND THE 
ONE-DIMENSIONAL CASE 

This Section lies somewhat outside of the main line of 



our development; the latter continues in Sec. VI. Here we 



consider our model in the vicinity of the NL deep in the 
low-temperature phase. In terms of the superspin chain 



with the Hamiltonian (4.7), in this phase we have /i <C fc. 
Then in the zeroth approximation we may neglect the 
h couplings completely. Then the chain (4.7) is broken 
into disconnected pairs of superspins. The Hamiltonian 
for one such pair is 



Hk = 4str AGAG = 4str AGAG, 



(5.1) 



where the coupling constant k (overall energy scale) was 
taken to be equal to 4 for later convenience, and we 
used the cyclic property of the supertrace. We can try 
to solve this Hamiltonian and hope to infer some infor- 
mation about the low temperature phase of our original 



model. However, we wish to sound a note of caution: 
we are considering a certain double limit of the original 
lattice Ising model, first the anisotropic limit, then the 
"low T" limit, h <^ k. It is not entirely clear that this 
really represents the low T limit of the nearly isotropic 
Ising model, where we pass to low T close to the NL, and 
perhaps then go to the anisotropic limit. 

We will make use of the realization of the representa- 
tions R and R in Fock spaces of unconstrained fermions 
and bosons. For simplicity we will work out the details 
for osp(3|2) only. In this case the necessary construction 
of R and R and the invariant products of superspins is 
given in Appendix H. From it we obtain 



with 



Hu = \J- r 



J = pp + ff + bW + bb. 



(5.2) 



(5.3) 



We anticipate that the eigenstates of the Hamiltonian 
Hk may have arbitrarily large bosonic occupation num- 
bers, and we may encounter convergence problems typi- 
cal in such cases. These are avoided , how ever, if we re- 
member the Lo term, discussed in Sec. IV B. As explained 
there, it plays the role of a symmetry-breaking regulator 
that picks a direction for ordering, similar to SI — 51 for 
the problem of two antiferromagnetically coupled su(2) 
spins. Thus, we add to our Hamiltonian the term 



The resulting Hamiltonian 

H = Hk + 



(5.4) 



(5.5) 



is identical to the one studied by Balents and Fisher in 
a one-dimensional localization problem (see Eq. (3.31) in 
Ref. ^8|). This is the problem of spinless fermions on a 
ID lattice with random hopping amplitudes described by 
the Hamiltonian 



(5.6) 



The continuum limit of this model gives left and right 
moving spinless Dirac fermions with random mixing be- 
tween them: 



(5J) 



where ^'(a;) is a two-component spinor field, are Pauli 
matrices, and the random potential is Gaussian with non- 
zero mean and variance: 



[{V{x) - Vo) {V{x') - Vo)] = 2DS{x - 



(5.8) 
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The generating functional for the Green's functions of 
this Hamiltonian at a given energy e + ir] may be su- 
persymmetrized in the standard way. After disorder av- 
eraging the X coordinate may be interpreted as imagi- 
nary time, and the two components of the fermion can 
be viewed as labehng two sites, which correspond to our 
spht sites. This leads to an effective quantum Hamilto- 
nian, which is exactly given by Eq. (p_^) with 



A = Vo/D, 



(5.9) 



The ID model with the Hamiltonian (5/7), and related 
models, have a long history, and most of the relevant work 
is concisely summarized in Ref. In particular, the 
density of states, for this problem was found for A = Orin 
1953 by DysoiB, and for arbitrary A by many authorsEa. 
The mathematically equivalent problem of diffusion i: 
ID random medium was studied by Bouchaud et a/I 
The density of states p{e) behaves at small energies as 



p(e) oc 
p(e) cx e 



1 



lln^l 

A-l 



A = 0, 



A > 0. 



(5.10) 
(5.11) 



In the superspin language the density of states p{e) is 
related to the expectation value of sorac. operator in the 
ground state of the Hamiltonian (p.5|)cJ. Namely, it is 
proportional to the staggered component /i2 — ^2 of the 
superspin (see Appendix |^) 



p{lu) cx (1 — n 



/ 



(5.12) 



This quantity measures the amount of the symmetry 
breaking in the ground state of two superspins. From 
the last equation it follows that the symmetry is spon- 
taneously broken on the NL (which is a point in the ID 
model) and below it. Moreover, below the NL, where 
A < 1, the density of states (the order parameter of the 
spin chain) diverges as w ^ 0. On the NL it is constant, 
and above the NL it vanishes as a power of uj. 

Because the SUSY representations are the same, we 
have in fact shown that in the ID off-diagonal disorder 
problem, there is a larger SUSY osp(27i -|- 1 | 2n) at the 
point A = 1. This has not been noticed previously to our 
knowledge. This suggests that such Nishimori points, 
lines, etc, may be common in some classes of random 
fermion problems. We also note here that in the ID 
classical RBIM, which of course has no finite T phase 
transition, there is a Ni shim ori point at which the cor- 
relation identities Eqs. (1.4) hold. That problem can be 
represented using fermions on one unsplit site with the 
Tii transfer matrices only, which are of the /i-coupling 
type, in contrast to the model considered here, and is 
easily solved in this language. 



VI. FINAL REPRESENTATION AND THE 
GENERALIZED MODEL 

In this Section we continue the general consideration 
of the RBIM problem. Here we focus our attention on 
the NL, that i s, w e consider the osp(2n-|- 1 1 2n)-invariant 
Hamiltonian ( |4.9| ). First we analyze and solve the prob- 
lem of finding a way to describe the term Npi in the 
spaces R, R on the split sites. The problem is solved 
by using another representation in a space JF", and the 
spaces can be viewed as representations of a larger SUSY 
algebra, osp(2n -t- 2 | 2n). Using only terms of the form 
of the two couplings we have already seen, we then in- 
troduce a more general nearest-neighbor superspin chain, 
and discuss its phase diagram, for reference in the follow- 
ing Sections. 

First let us note that Eq. (^^) is somewhat schematic. 
Let us again consider the fcrmionic replica formalism 
with n ^ instead of SUSY. The term N'g^ is then re- 
placed by N'p-. According to our general discussion of 
how to map operators in into !F' (see Sec. Ill), 



a 



(6.1) 



is correct as it stands in T' . Even though the operator 
Np^ is perfectly legitimate, it does not admit any simple 
expression in terms of the so(2n -I- 1) generators Gi and 
Gi. We would like to write it as a sum of products of 
operators in the spinor representations R, R on the split 
sites. Of course, individual fermion operators ^ai, Vai do 
not commute with the constraint, and cannot be used. 
Instead they must be replaced by Z2-invariant operators. 
As explained in Sec. Ill , we can find operators in J-' with 
the anticommutation properties of the fermions in T for 
the components other than the zer oth. These are fermion 
bilinears times a string; see Eq. (3.9). A general proof 
that it is impossible to find a set of operators with the 
anticommutation relations of the full set of real fermion 
operators rjai and S^ai in the space is to notice that they 
should form a Clifford algebra with 2N{2n + 1) genera- 
tors, where N is the number of unsplit sites in the chain. 
This Clifford algebra has a single non-trivial representa- 
tion of dimension 2^^^"+^). This space is the same as 
an unconstrained Fock space for 2n -I- 1 complex fermion 
operators at each unsplit site, i.e. J^' but without the 
constraints. The total number of states in JF' is only 
2^^", because of the N constraints. So the operators 
we require cannot have the anticommutation relations of 
free fermions for all the sites. Indeed, in our grading on 
single fermion operators are even, and so would be 
expected to obey commutation relations from a SUSY 
point of view. In J-\ there are corresponding fermion 
bilinears, like Eq. ( ^.9[ ) but without strings, and these do 
commute on different sites (as mentioned already in Sec. 
11). 

We can also try the unconstrained representation. 
Then again, we can represent each of the 2n real fermions 
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rjai and ^ai on each split site using Eqs. ( 4.1C| ), and the 
resulting Clifford algebra for 2N split sites yields the cor- 
rect number of states. This description carries over easily 
to the SUSY version. But the above proof shows that no 
matter what strings or other factors we introduce into a 
construction of operators, we cannot produce the anti- 
commutation relations for 2n -\- I real fermions at each 
split site, and we are no nearer writing N'g^ as a product 
of simple expression in R and R. What we would have to 
do is map the problematic part ngi of the operator back 
from to F . Because of the constraint in the former 
space, the resulting operator (still in the fermionic replica 
formalism) must equal 1 when Npi is odd, when Npi 
is even, or similarly for Ngi in the SUSY formalism. It 
is not clear how we would write this as a coupling of the 
two split sites at i. 

There is nonetheless a way out of this problem, mo- 
tivated by the following observation. If we consider a 
single spinor representation R of so(2n -t- 1) (thus, in the 
fermionic replica formalism once more) , then it is in fact 
possible to find operators with the anticommutation re- 
lations of the real fermions ^q. These are the generators 
of a Clifford algebra with an odd number 2n -I- 1 of gener- 
ators, which has an irreducible representation of dimen- 
sion 2" (the familiar 2x2 Pauli matrices are the case 
n = 1). The commutators W^a,^t\ of these operators are 
the generators of so(27t. -|- 1), as we have already seen. 
The operators <^a transform as a vector of so(2n + 1). 
If we now consider these so(27i + 1) generators together 
with the fa (divided by -^72), then we may use the fact 
that the commutators of these operators (or matrices) to- 
gether obey the relations of the generators of so(2ri, -I- 2), 
and the spinor R can be identified with one of the two 
distinct irreducible spinor representations of dimension 
2" of so(2n -I- 2). This construction can also be applied 
to the representations R [which for so(2n + 1), though 
not for osp(2n + 1 1 2n), is isomorphic to R]. This con- 
struction also extends easily to the many-site problem, 
by taking the operators now replacing fa to commute on 
different sites. We may therefore write down our term 
Np^ as a sum of products of bilinears of these operators, 
and this can also be extended to the SUSY construction, 
using operators with the relations of g^, r„, on each site, 
but which anticommute on different sites. 

Thus, we have learned that our spaces of states i?, R at 
alternate sites can be viewed as irreducible spinor repre- 
sentations of osp(2n -t- 2 I 2n) (but note that this algebra 
is not a symmetry of our Hamiltonian or transfer ma- 
trices so far). There are two inequivalent lowest-weight 
spinor representations Re, Ro (for "even" and "odd") of 
osp(2n + 2 I 2n), in which all states can be assigned posi- 
tive norm-squares. We can identify R with, say. Re- Sim- 
ilarly, R can be identified with a highest-weight spinor 
i?e, which is dual to Re, and in which the inner product 
is indefinite, since as we have seen states with an odd 
number of b bosons have negative squared norms. 

Viewing the spaces of states in this way, we can give yet 



another explicit construction, with which we can finally 
write the operators N'p^ and N'g^ in a simple way. It 
is convenient to keep much of the notation the same as 
before. We introduce additional complex fermions /oi, 
/oi to the set f^t of Eq. (4.1C), and define a space T" — 
®iT'l consisting of the states on the split sites with an 
even number of fermions plus bosons: 



^fi + flJoi = even, 
"/i + flJoi = even. 



(6.2) 
(6.3) 



It is clear that such states are in one-one correspondence 
with those in the unconstrained representation. The con- 
struction of the correspondence is similar to that for the 
states in the spaces !F and f in Sec. III. All the states 
can be obtained from the vacuum, which is the lowest- 
(highest-) weight state in Re (Re), by the action of the 
bilinears in the creation operators. Then in addition to 
Eq. ( pTO| ) we also define 



/o 



V2 



foi — 



V'oi 



V2 



(6.4) 



Now the so(27i-|-2) generators on a single site are replaced 

by 



2 Kaij fbi] I 



2 Koii ^ai], 



(6.5) 



where the first set, again, spans the subalgebra so(27i+l), 
and the second set transforms as a vector under this sub- 
algebra. There are similar expressions for Re- We em- 
phasize that the operators f a i , Coi i Vai , v'oi obey canonical 
anticommutation relations, while qai, r^i obey canoni- 
cal co m muta tion relations, of the same form as in Eqs. 
( p. 17 ), (2.24), (with negative norm states appearing in 
connection with the r's, see Eqs. ( 4.11 ) and following). 
Our choice of grading is again equivalent in the con- 
strained subspaces J^-' to their natural grading as sub- 
spaces of Fock spaces. Finally, in the representation in 
T" , the operator Ng.^ undergoes the replacement 



(6.6) 



These results may also be established by passing directly 
from the (averaged) unconstrained representation to the 
final representation, by using a substitution similar to 
Eq. (3^), but applied here to the split sites. 

We can now organize the g enerator s of osp(2n + 2 | 2n) 
in superspins, similar to Eqs. (O, O), and including the 
additional odd generators: 



G' = 




Co6 

^a£,b — 2^ah 



VaVb 



Co 9/3 

fag/3 ^ I, (6.7) 

iJa^q^qfi — 2^af3 , 

Vari3 . (6.8) 
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Note that these osp(2n + 2 | 2n) superspins contain the 
original osp(2n + 1 1 2n) superspins G and G as subma- 
trices. The odd generators are those containing an odd 
number of fermion operator factors, or equivalently an 
odd number of boson operator factors. With the help of 
the osp (2n + 2 | 2n) superspins, both terms in the Hamil- 
tonian (4.9) may be written in a unified way: 



hN' = lim ( str CG'CG' + - 



(6.9) 



kstiG^Gi+i ^ li-m ( str CG-CG-+1 + 2 ) ' (^-lO) 



where we have introduced a An - 
matrices of coupling constants 



' 2-dimensional diagonal 



G = C{h, k) = diag(/ifc-i/2^ A:i/2lL4„+i), (6.11) 

and the supertrace str in this space is defined in the same 
way as the previous str . These two terms represent two 
different osp(2rt-|-l | 2n)-invariant products of two subsets 
of the osp(2n + 2 | 2n)generators. It should be clear that 
the representation in T" can also be used off the NL, by 
giving certain terms different coefficients. 

It is now natural to consider a generalized Hamiltonian 

Y,isii-C{hA,kA)G',C{hA,kA)G', 



+ stYC{hB,kB)G',C{hB,kB)G',+,] 
h-A + h-B 



(6.12) 



parametrized by four coupling constants. In such a 
Hamiltonian both types of osp(2n+ 1 1 2n)-invariant cou- 
plings appear on every bond between the split sites. In 
addition, they are staggered between the two sublattices 
of bonds A and B of our chain. Our NL Hamiltonian 



(4.9) is a particular extreme limit, where on alternate 
bonds on e or t he other coupling is zero. It is obtained 
from Eq. ( 6.12| ) for the special values of the parameters 



llA = /i, kA^ 0, 



hu =0, ku = k. 



(6.13) 



We believe it may be helpful to consider these more gen- 
eral models, since they are so closely related to that for 
the RBIM, and use only couplings that appear anyway in 
the RBIM case; however, we emphasize that it may not 
he possible to obtain these models as anisotropic limits of 
random fermion or network models. When 



hA = kA, 



hs = kf 



(6.14) 



the model is invariant under the whole of osp(2n-|-2 | 2n). 
We should note that in principle we can also consider this 
generalization in the discrete imaginary time model, and 
also off the NL, where however the breaking of the sym- 
metries would lead to twice as many parameters. The ad- 
ditional parameters would generalize A, A* in Sec. |IV A , 
and there would be one for each of kA, ks, Iia, hB, a 



A:. = 00 




FIG. 3. A possible phase diagram of the generalized Hamil- 



tonian (6.12), discussed in detail in the text. 



total of eight parameters in the HamiltonLajji. In particu- 
lar, another model due to Cho and Fisherll3 fits into this 



general description, as we will see in Sec. VII. 

The Hamiltonian ( 6.12| ) contains four parameters, but 
since the overall energy scale is unimportant, the phase 
diagram can be plotted in terms of the three indepen- 
dent ratios of parameters. We will consider only positive 
values of all couplings, though negative values may also 
give well-defined models. The phase diagram can then be 
drawn in a symmetrical manner in a three-dimensional 
tetrahedron, as a portion of projective space (see Fig. |^). 
Each face of the tetrahedron is defined by one of the four 
parameters vanishing. The opposite vertex is where that 
parameter goes to infinity, or equivalently the other three 
all go to zero. 

The edges of the tetrahedron correspond to models 
with two vanishing couplings. For example, the verti- 
cal edge connecting the ver tice s Ha — 00 and ks = 00 
represents the Hamiltonian ( O ) for the RBIM on the NL 
(to avoid confusion, recall that the whole discussion is a 
generalization of the NL, since all the models have the 
larger osp(2n-|- 1 1 2n) SUSY). There is another such line 
represented by the horizontal edge connecting the ver- 
tices hB = 00 and kA — 00. These two Hamiltonians are 
related by a reflection through a lattice site. Such an op- 
eration is thus a symmetry of the whole diagram, which 
interchanges A with B. The line kA = kB, hA = hB is 
invariant under this operation, and the operation acts as 
a 180° rotation about this line. On each NL, there is a 
multicritical point, N and its image N' . 

The edges where both non-zero couplings are on the 
same sublattice of bonds (e.g. A) represent the two ex- 
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treme cases of fully-dimerized chains, which have a gap 
in their energy spectrum. By analogy with other anti- 
ferromagnetic (super-)spin chain models, we expect that 
the regions adjacent to these lines are also gapped phases. 
There must be at least one phase transition between these 
two extremes. One way to see this is to consider a chain 
with open ends, and an even number of split sites. In 
one phase the dimers extend all the way to the ends of 
the chain, in the other a single superspin is left unbonded 
with a neighor at each end. This corresponds to a chi- 
ral edge degree of freedom in the 2D lattice model. A 
phase transition must occur to change the number of such 
boundary spins or edge channels, assuming these survive 
off the edges of the tetrahedron. We will assume that, 
as expected on the NL in the RBIM, there is a single 
transition between the two phases. Then there must be 
a phase-boundary surface between those two edges, indi- 
cated schematically (since its exact position is unknown) 
by the shaded surface in Fig. ||. The points N and N' are 
two vertices of this rhomboidal surface, which also con- 
tains the line of reflection symmetry. However, we note 
that an intermediate phase, in place of some portion of 
the critical surface, is also possible, though this is not 
expected on the NL in the RBIM. 

The two other edges of the tetrahedron are where ei- 
ther only the fc's are nonzero, or only the /I's, and we 
denote these models "fc-only" and "ft,-only" . They inter- 
sect the phase boundary (if there is a unique transition on 
these edges) at points labeled K and H in Fig. ^ (no con- 
fusion should result from this notation). In these models, 
the reflection symmetry and the assumption of a single 
transition implies that K is kA — ks, and H is Ha — hs 
(and other parameters zero). 

The tetrahedral phase diagram also contains the line 
given by Eq. ( |6.14 ) where the generalized Hamiltonian 
(3.12) has the osp(2n-|-2 | 2n) symmetry. This line, shown 
dotted in Fig. ^, intersects the critical surface at a critical 
point O (black dot), where all four couplings are equal. 
Again, this point is unique if we assume there is a single 
transition on this line; it is kA — ks ~ hA ~ hs- 



VII. CHO-FISHER, fc-ONLY, AND /i-ONLY 
MODELS 

In this Section we consider a model studied numerically 



by Cho and Fisher in Ref. 13, This is a network model, 
similar to the Chalker-Coddington networlalil describing 
the integer quantum Hall transition, but with only real 
matrices, and was intended to represent the RBIM prob- 
lem. We show that the Ising model can be represented 
exactly as a network model, and that the Cho-Fisher 
model does not represent the RBIM. Instead, it can be 
mapped to some of the generalized Hamiltonians without 
enhanced SUSY, introduced in Sec. VI. 




The Cho-Fisher model is a network model, intended 
to capture the universal aspects of the point N in the 



FIG. 4. Relation of the Ising model and the network model. 
Ising spins are located at the open circles, and bonds are 
shown dotted. Solid lines with arrows form the "medial 
graph", on which the network model is defined. Examples 
of nodes on each of the two sublattices, corresponding to the 
horizontal and vertical bonds, are labeled A, B, respectively. 
In this paper, we consider periodic boundary conditions in the 
vertical direction, and free in the horizontal, as shown here. 

RBIM, which can be viewed as a generalization of the 
ID model discussed in Sec. 0. It was constructed as a 
generalization of the model whose action is given in Eq. 
( ^.7[ ) , in which the two components of the fermion are re- 
placed by any number of sites in a ID chain, with random 
nearest-neighbor hopping that generalizes the term in 
Eq. ( |5.7| ), and, in general, different parameter values Vq 
and D on alternate bonds in the chain. Then use of repli- 
cas or SUSY to perform the disorder average leads to a 
generalization of th e quantum (super-)spin Hamiltonian 
Hk + H^, Eq. (5^) in Sec. (we will disregard here the 
regularising term H^), in which the same form of cou- 
pling appears for each pair of nearest neighbors, but with 
the coefficient of Hk taking two values k, k' on alternate 
bonds, and similarly for A, A'. We emphasize that at 
this stage we are using the unconstrained representation 
of the space of states of the chain. Cho and Fisher spe- 
cialized to the case k — k' (i.e. node independent disor- 
der strength), and went back from their time-continuum 
model to a discrete-time (network) model, similar to that 
in Ref. |l^, in order to perform numerical calculations. 
They claimed that their network model has a multicrit- 
ical point in its phase diagram with critical properties 
remarkably similar to those of the multicritical point on 
the NL. In particular, the critical exponents along the two 
scaling axes near the multicritical point were found to be 
fairly close numerically to the ones known for the point 
N in the RBIM pmm the high-temperature expansion of 
Singh and Adlerll3. Also, simulation in Cho's thesial3 of 
a network model that corresponds precisely to the RBIM 
(as we will explain) gave similar values. 

Now, we wish to point out that it is possible to relate 
network models more directly to the transfer-matrix for- 
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mulation for fermions, as in the Ising model. First we 
describe the network models. A portion of the network 
is shown in Fig. ^, where the solid lines with arrows are 
where the particles propagate. The particles propagate 
in discrete time, at each time step moving to the next 
link in the "forward" direction shown by the arrows, and 
therefore turning either left or right at each time. The 
evolution is described by a unitary S-matrix which gives 
the amplitudes for turning either right or left at each 
nodetil. This can be replaced by a one-particle trans- 
fer matrix, which adds one row of nodes to the system, 
evolving the wavefunction of the particle upwards in the 
figure. In the Cho-Fisher model, the one-particle transfer 
matrix for one node has the form 



M 



cosh 9 sinh 9 
sinh 9 cosh 9 



(7.1) 



The parameter 9 for each node is random, taking val- 
ues ±1^1 with probabilities 1 — p, p, independently. Also, 
the magnitudes of 9 can be staggered, taking different 
values \9a\^ \0b\ on the two sublattices of nodes labeled 
A, B in Fig. ^. The so-called isotropic case is where 
sinh \ 9a\ sinh \ 9b\ = 1. This leaves a one-parpiieter fam- 
ily of models; in the original network modelij, a transi- 
tion occurred when \9a\ = the "self-dual" point. 

To exhibit a relation with the Ising model transfer ma- 
trices, we use a second-quantized formulation of the net- 
work, as a noninteracting fermion field theory. The evo- 
lution in the imaginary-time (vertical) direction is de- 
scribed by a transfer matrix constructed from the one- 
partijeie one. We can write this by drawing on earlier 
workO. Though the latter was on a different model, the 
basic Eq. (2) in that work is applicable for any transfer 
matrix, with matrix elements 



a P 
7 5 



(7.2) 



Thus in our case, a — 5 ^ cosh^, = 7 = sinhf?. Using 
only one species of fermions, and dropping the bosons in 
Eq. (2) of Ref. ^ since we will not be averaging here, we 
replace /i in Ref. |^by /, /2 by /', (where / and /' obey 
canonical anticommutation relations) and obtain 

V exp[tanh6l(/t/' + ft/)] : (cosh6')"^+"/' . (7.3) 

Here the colons : . . . : indicate normal ordering with de- 
struction operators to the right. Then, after making the 
particle-hole transformation p = /', we can prove the 
identity 



(7.4) 



by verifying that all matrix elements of the two ex- 
pressions are equal. But this now has the form of the 
fermionic representation of the squared Ising model (i.e., 
n = 1), as in Eqs. (2.32), (4.13) (dropping the bosons 
in the latter), up to constant factors in the vertical case. 



This means that the split sites on a single row correspond 
to one row of links of the network. The relation of the 
original Ising lattice and the network model is as shown 
in Fig. H in particular, the two sublattices of nodes A, B, 
correspond to horizontal and vertical bonds respectively. 
The relationship of second-quantized transfer matrices 
holds true for arbitrary values of 9 at each node, and 
also remains true when bosonic partners are introduced 
in preparation for averaging. 

The important corollary to this is that for the transfer 
matrices of the sort appropriate Jar the horizontal bonds 
(labeled A in Fig. |), we haveN 2K = 9. The Cho- 
Fisher network model takes the parameter 9 at the nodes 
to have independent random signs. Hence i t is pr ecisel y 
equivalent to the use of transfer matrices (2.9), ( 2.10 ), 
with the binary distribution of the type ( |1.!^ ) for both 
the horizontal couplings Ki^i^x and the dual Ki^i^y to 
the vertical couplings. Note that the Cho-Fisher model 
is not in fact isotropic, even when the (magnitudes of 
the real parts of the) K's, and the probabilities p, on 
the horizontal and vertical links are the same, which is 
what we termed isotropic above; this is because of the 
way the random signs are introduced. Another popular 
parametrizatian for the network models uses the S-matrix 
at each nodeliJ, where the S-matrix is a real orthogonal 
matrix in the present case, with one of its off-diagonal 
matrix elements (say, the amplitude for turning right) 
denoted t ~ simp. In this case the equivalence to the 
Ising model squared is tanh 2K = sin (/) for the horizontal 
couplings. 

Since negative dual couplings K correspond to complex 
Ising couplings K, the Cho-Fisher model does not faith- 
fully reproduce the RBIM with ±K couplings. Instead, 
in replicated fermion language, both types of bond are 
represented after averaging by transfer matrices of the 
horizontal type. One might imagine that this is the fc- 
only model, with parameters Xa, Ab included, so that 
the osp(2n + 1 1 2n) SUSY is present when A = A' = 1. 
But in fact, carefully following our mapping leads to dif- 
ferent forms for the two bonds. A, B. It is necessary 
once again to use the final {J-") representation, and to 
pass to it directly from the unconstrained representation. 
We find that, while the bonds corresponding to the hori- 
zontal bonds in the Ising model involve the real fermions 
Coij Voii in the A:-coupling terms as in Eq. ( 6.10 ), for the 
vertical bonds those fermion operators are replaced by 
Coii v'oi- For A or A' = 1, these terms are invariant 
under an osp(2n -|- 1 | 2n) SUSY, but these are distinct 
osp(2n + 1 1 2n) sub-superalgebras of osp(2n -I- 2 | 2n) in 
the two cases, and so the Hamiltonian does not possess 
a global osp(2n + 1 \ 2n) SUSY (though there is of course 
still osp(2?i I 2n)). These models therefore lie elsewhere 
in our space of fully-generalized Hamiltonians with (in 
general) only osp(2n | 2n) SUSY. 

We note here that by rotating the Cho-Fisher network 
by 90° , we obtain after averaging (using the !F" represen- 
tation) a model which resembles the /i-only model, but 
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FIG. 5. The critical surface of the generalized Hamiltonian, 
after reduction using symmetry. The arrows indicate our sug- 
gested RG flows to a fixed-point theory with the same larger 
SUSY as the model at O. 



again has different osp(2n -I- 1 1 2n) SUSYs for the two 
types of bond. Because taking the anisotropic hmit usu- 
ally does not affect the universality class of the critical 
phenomena, the resulting spin chain model should have 
the same critical phenomena as the one described above 
for the Cho-Fishcr model. 



VIII. FIXED POINTS AND NON-LINEAR 
SIGMA MODELS 

In this Section, we first speculate that a single fixed 
point, or universality class controls much of the phase 
boundary in the tetrahedral phase diagram. Then we 
discuss the nonlinear sigma models that are related to 
our spin chains, and calculate, at weak coupling, the RG 
beta functions for the coupling constants. The results 
support the hypothesis of a flow towards the higher SUSY 
as at point O in the phase diagram. Finally, we discuss 
the nonlinear sigma (and related) models for the more 
general, lower SUSY (osp(2ri | 2n)), or class D, random 
fermion problems. 

For the generalized Hamiltonian with osp(2n + 1 1 2n) 
SUSY, we argued (on the assumption that there is a 
single transition surface) that the phase boundary is a 
rhombus, and further there is a reflection symmetry in 
the superspin chain, which, on this surface, acts as a re- 
flection. Then the phase boundary is a triangle, with 
points N, H at the vertices, as shown in Fig. ^ The 
point O, at which the model has the larger SUSY algebra 
osp(2n + 2 I 2n), is now at the middle of one side of the 
triangle. 

Under the RG, the osp(2n -I- 2 | 2ri)-invariant model 
must flow to a critical quantum field theory (presumably, 
a conformal field theory) which also has the larger SUSY. 
Other models, represented by other points in Fig. |^, such 
as iV, ivT, H may flow to some other fixed point theories 



of lower (osp(2ri + 1 1 2n)) SUSY, and it is of course the 
fate of N that concerns us in the RBIM model problem. 
If we must make a guess as to the structure of the flows 
and flxed points, the simplest guess is the one that in- 
volves the fewest flxcd points. Since there must be a 
fixed point theory with osp(2n + 2 | 2n) SUSY, the sim- 
plest guess is then that the whole critical surface flows to 
this fixed point. This is schematically illustrated in Fig. 
I by the arrows, which are intended to indicate that all 
models fiow to the fixed point corresponding to O (note 
that the models at N , K, H, and O are not themselves 
fixed points of the RG). If correct, this would imply that 
the critical exponents are the same at all points in the 
critical surface shown in Figs. ^ or ^. In particular, 
N, K, and H would have the same exponents. We can 
also imagine other scenarios in which iV, K and H flow 
to a common flxed point, or to different ones, that do 
not have osp(2n -I- 2 | 2n) SUSY. It is certainly possible 
that (one or both) perturbations away from O on the 
surface shown are relevant; we are suggesting that they 
are both irrelevant. In the absence of any understanding 
of the conformal field theory of the osp(2n -I- 2 | 2n) or 
other fixed points in this system, we cannot prove or dis- 
prove our suggestion. There are, however, other systems 
in which an analogous effect occurs, as we will discuss 
below. 

Now we introduce the nonlinear sigma models that 
should correspond to the superspin chains, and should 
have transitions in the same universality class or classes. 
First we utilize a standard relation between antiferro- 
magnetic spin chains and nonlinear sigma models (see 
Refs. ^ for a fairly general discussion). We define an 
antifcrromagnetic (super-)spin chain as having an irre- 
ducible representation at each site, alternating between 
some (say) lowest-weight representation R and its dual, 
say R. The Hamiltonian should be something close to the 
Heisenberg form which is the invariant bilinear form in 
the generators of the symmetry algebra, with the antifcr- 
romagnetic coupling that for a single pair of spins leads to 
the singlet ground state (possible because we chose the 
dual representations). Then the correspondence states 
that there is a nonlinear sigma model with a certain tar- 
get (super-)manifold, which can be obtained from the 
representation R. The manifold is the same coset space 
that appears as the coadjoint orbit of R, or in coherent- 
state path integral constructions of R. Put simply, this 
is the manifold swept out by acting on either a lowest- 
or highest-weight of R with all possible (super-)group 
elements. The long-wavelength action of the nonlinear 
sigma model in 1 -I- 1 dimensions contains only terms 
allowed by symmetry with two derivatives. These com- 
prise the usual "kinetic" type terms, and also possible 
"6>-t erms " (this is not the same parameter 9 we used in 
Sec. VH). The derivation is controlled by considering a 
sequence of representations R with the lowest weight go- 
ing to infinity in the weight space, like the size of the 
spin in SU(2) going to infinity. Then the reciprocals of 
the kinetic couplings have magnitude proportional to the 
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lowest weight, so the nonhnear sigma model is weakly 
coupled and meaningful in the semiclassical limit. Also, 
in the absence of staggering of the couplings in the spin 
chain, 9 is proportional to the lowest weight with a coef- 
ficient TT in a suitable normalization (which is such that 
the bulk physics is periodic when 6^6 + 2tt) . A 9 term 
exists and is nontrivial whenever the second homotopy 
group 7r2 of the target manifold is nontrivial. More gen- 
erally, a 0-term involves a two-form on the target man- 
ifold (i.e. a magnetic field for a charged particle moving 
on the manifold) that is invariant under the symmetry, 
and this always exists in this construction because it is 
part of the coherent-state construction of the represen- 
tation R also. Noncompact factors in the manifold are 
topologically trivial, but the term described always pro- 
duceSphpundary effects related to a boundary spin or edge 
statesE3. The nonlinear sigma model that results from 
this correspondence in many cases has a phase transition 
at = TT, when the target manifold has nontrivial 7:2- 

In our case, the target manifold for our general 
models on the NL would be, for fermionic replicas, 
S0(2n + 1)/U(n), or in the SUSY formalism, 0Sp(2n + 
l|2n) /U(n|n). The precise meanings of these coset spaces 
should be defined as the orbits of our spinors R. The 
group in the denominator arises in each case as the invari- 
ance group of the lowest weight state. In the SUSY case, 
U(n|7i) corresponds to the superalgebra gl(7i|n) which 
leaves the vacuum at each split site invariant, and the 
notation indicates that the ordinary group it contains 
is the compact form, U(n)x U(n). The manifold un- 
derlying the supermanifold is thus S0(2n -f 1)/U(n) x 
Sp(2ri, R)/U(n), where the latter factor is non-compact. 
These (super-) spaces are homogeneous spaces, but not 
symmetric (super-)spaces. This implies that, in a gen- 
eral nonlinear sigma model for these target manifolds, 
there is more than one coupling in the kinetic termsc2l — 
in fact, there are two (see Appendix]^). The kinetic term 
is constructed from the metric on the target manifold, 
and this metric is not unique, up to a constant factor, 
unless the manifold is a symmetric space. Otherwise it 
is a sum of two or more pieces. However, as manifolds 
(without a choice of metric), these spaces are the same as 



2n)/\]{n+l\n)] 
ni) . This is con- 



S0(2n + 2)/U(n + 1) [resp. 0Sp(2n + 2 : 
(for the former, this is discussed in Ref. 
nected with the fact that, as graded vector spaces, R is 
the same as the representation R^. The latter manifolds 
are symmetric (super-) spaces, and there is a unique ki- 
netic coupling. The two couplings in the former point of 
view are related to the k and h terms. When k and h on 
each nearest-neighbor bond are equal, the higher symme- 
try implies that the two coupling constants in the kinetic 
terms are such that the osp(2n -f 2 | 2n) invariant kinetic 
term is obtained. At the same time, 9 can be varied by 
staggering the couplings. Counting parameters, there are 
four in the spin chain as we have seen, one of which is 
the overall energy scale which can be ignored, and also 
the magnitude of the lowest weight. The latter controls 
the magnitude of the kinetic coupling constant. Ignoring 



that, the nonlinear sigma model has in general one ratio 
of kinetic couplings, and 9, one parameter less than in 
the spin chain. However, so far in the nonlinear sigma 
model wc have assumed Lorentz invariance, whereas in 
the continuum limit of the spin chain this need not hold. 
In fact, in the sigma model, the same symmetry consid- 
erations imply that each of the two pieces of the kinetic 
term can be further divided into two terms, which are the 
second derivatives in the two orthogonal directions in the 
2D space, and if Lorentz invariance is not required, all the 
terms can have different coefficients. That is, a different 
velocity can. apply to the two kinds of spin waves in the 
spin chainH; the two kinds of spin waves correspond to 
the decomposition of the small fluctuations around the 
perturbative vacuum into two irreducible representations 
of the symmetry group U(n) [or U(ri|n)] (see Appendix 
^). Since the overall energy scale, or one velocity, is a re- 
dundant parameter, this leaves one additional parameter, 
the ratio of velocities, as required. Note, however, that if 
we instead consider the continuum limit of an isotropic 
network model, such as the RBIM, only one velocity can 
occur. Therefore we would expect that the universality 
classes of the transitions would be isotropic and have a 
unique velocity. 

Next, it is natural to raise the question of the RG flow 
of the two couplings in the osp(2n-|- 1 1 2n)-invariant non- 
linear sigma models. At least in perturbation theory, the 
corresponding beta functions that describe such a flow, 
are the same as those for the models with ordinary target 
manifolds S0(2n-|- 1)/U(n), with n — > 0. In Appendix |Q, 
we have computed the beta functions to one-loop order, 
neglecting the possibility of more than one velocity (in 
some analogous situations, thej'atio of velocities has been 
shown to renormalize to onecHI, and we expect the same 
to occur here, as also argued in the last paragraph) . The 
metric on the target manifold is parametrized by two pa- 
rameters, say 7?!, 772, which appear linearly in the metric, 
and the inverses of these parameters are the couplings 
which are small at weak coupling, where perturbation 
theory is valid. The perturbation expansion has l/?7i, 
l/?72 for each propagator, 771 or 772 for each interaction 
vertex, so naively one ends up with a perturbation ex- 
pansion in, say, l/?7i, with each term containing powers, 
positive or negative, of the ratio of 771 to 772. In fact, only 
non-negative powers of this ratio appear in the beta func- 
tions; see App. p. The net power (counted with signs) 
of 7;j^^'s and 77^ 's corresponds to the number of loops 
in the Feynman diagrams, as usual. For convenience we 
will use the parameters 
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2xg' 



(8.1) 



In this parametrization, .t = 1 is the point with osp(27i -|- 
2 I 2n) SUSY (or so{2n + 2) in the replica version). Be- 
cause of the higher symmetry at a; = 1, that line should 
flow onto itself under RG. Then our one-loop result for 
the RG flows is (/ is the logarithm of the length scale, as 
usual) 
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FIG. 6. Sketch of perturbative renormalization group flows 
for the couphngs g, x, including two-loop effects, for the non- 
linear sigma model with target space OSp(2n-|- l|2n)/U(n|n), 
including a typical flow line for g ^ 0, x 0, 1. 



(8.2) 
(8.3) 



At a; = 1, the one- loop result is zero, so this line is a line 
of fixed points, to this order. This agrees with the one- 
loop result for the osp(2n + 2 | 2n)-invariant model; the 
beta function to two-loop order, obtained as the ti — + 1 



limit of that for the S0(2n)/U(n) model (see Eq. (|C10|) ), 



IS 



0{g'). 



(8.4) 



Thus, it vanishes to one-loop order, but not to two-loop 
order. At two-loop order, g flows towards strong cou- 
pling (see Fig. ^). The one- loop flows for x 1 take 
X closer to 1, and (except for a; = 0) the flows starting 
at (7 ^ never reach g = 0. Instead they flow to the 
region a; ~ 1 where the one-loop terms vanish, and the 
two-loop term cannot be neglected. Since the two-loop 
term in dg/dl at x = 1 is positive, all flows from weak 
coupling eventually go towards large g, with x approach- 
ing 1, except when x = 0. On the latter line, g flows 
towards weak coupling. As discussed in App. |^ (but 
here in SUSY language), on this line the larger SUSY 
has spontaneously broken, and the system is described 
by the OSp(2n|2n)/U(n|n) nonlinear sigma model [there 
is an additional global degree of freedom, described by a 
point on 0Sp(2n + l|2n)/OSp(2n|2n), a "supersphere" , 
on which the larger SUSY algebra acts] . Flows that be- 
gin at small, nonzero x eventually go to strong coupling. 
This generates very large crossover lengths, due to the 
very slow flows near x = 1, g = 0, where the flrst nonzero 



term is at two-loop order; the length scale at which g be- 
comes of order one is of order exp[l/(2(7o2;o) + l/(85o^o)]' 
in units of a short distance cutoff, where go and xq are 
the bare values of g and x, and it was assumed that xq 
and goXo are both small (see App. ^). 

As usual, the perturbative results do not depend on 9, 
but such dependence can be expected nonperturbatively. 
For the sigma model with osp(2n -I- 2 | 2n) SUSY, and 
for X ^ 0, the flows go towards strong coupling, and it 
is highly plausible, based on our experience with transi- 
tions (such as the integer quantum Hall transition) with 
such behavior of the couplings, that there is a unique 
fixed point at strong coupling g and at 9 — tt (mod 
27r) . Hence we expect that the spin chain at point O has 
the same critical theory as the osp(2n -I- 2 | 2n)-invariant 
sigma model. It is also quite plausible, based on the be- 
havior of the flows, that at least models that map onto 
the nonlinear sigma models with osp(2n + 1 | 2n) SUSY 
also flow to the same critical theory, when 9 = 7: (mod 
27r). It is still possible that in our generalized spin chains, 
points other than O do not flow to the osp(2n -I- 2 | 2n) 
critical theory, but it is plausible that there is a nontriv- 
ial neighborhood of O on the critical surface that does. 
This possibility may seem more plausible if we point out 
that in some other cases (without SUSY), a similar phe- 
nomenon is believed to occurc3. It is of course less clear 
that distant points such as N, K, and H flow to the 
same theory. Since the spin chain models typically start 
at bare couplings of order 1, we can almost rule out any 
flow to the weak-coupling regime of the osp(2n-|-l | 2n) or 
osp(2n + 2 I 2n)-invariant nonlinear sigma models, (anal- 
ogous to that,in_the lower-SUSY osp(2n | 2n) nonhnear 
sigma modelEaEll) because that regime is not stable un- 
der the RG. A flow to weak coupling is only possible by 
tuning a parameter, corresponding to putting x = or 
g = in the weak-coupling analysis. A natural guess is 
that the /i-only models might satisfy one of these con- 
ditions. This might even occur for a range of values of 
the staggering, corresponding to changing 9 away from tt 
(mod 27r) in the nonlinear sigma model, since the value 
of 9 is irrelevant (or formally, exactly marginal) at weak 
coupling. We have been unable to see why any of these 
models should satisfy such a condition exactly. However, 
it may be that one of them lies close to x = 0. In that 
case, the RG flows take them close to g = 0, and since the 
flows to strong coupling pass near x = 1, g = 0, where the 
first nonzero term is at two loops, the crossover length 
could be very large. That is, very large systems would be 
needed to see the true asymptotic critical behavior. An- 
other possibility is that these models lie at bare values 
X > 1, g very small, which again yields a large crossover 
length. Since we have a two-parameter space of critical 
models, we may expect to be able to tune g or x small 
somewhe]f&|in this space. However, arguments presented 
elsewhereE3 show that the RBIM, and hence the point N, 
cannot flow to the weak coupling region. 

The leading alternative scenario seems to be that al- 
though most points in the phase diagram flow to the 
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osp(2n + 2 I 2n)-invariant fixed point, the multicritical 
point N on the true NL, may be a distinct universality 
class, and the perturbation off this point in the phase 
diagram may be a relevant one that causes a flow to 
the osp(2n + 2 | 2n)-invariant fixed point. Clearly, we 
cannot answer here the question of which of these sce- 
narios is correct. But the self-duality apparent in the 
osp(2n + 2 I 2n)-invariant model at its critical point, as 
manifested by the reflection symmetry about a split site 
in the chain, and the significant lack of it in the RBIM, 
which instead has the special riisymmetry" property that 
the Ising couplings are all reaE3, suggests that this alter- 
native scenario may be correct. 

There is one further point to make about the spin 
chains and nonlinear sigma models, that applies off the 
NL. In that case the SUSY is broken to osp(2n | 2n). It 
will be convenient here to make use of the replica for- 
malism, in which the language and notation are simpler 
and standard, but the ideas extend also to the super- 
groups, since the additional bosons and Sp(2n,R) sym- 
metry, and the odd generators, do not change the form 
of the argument. In the higher SUSY Hamiltonians, in- 
cluding that for the NL, the global symmetry group can 
be seen to be S0(2n -|- 1), and making any of the A pa- 
rameters ^ 1 reduces the symmetry group to 0(2n), not 
just S0(2ri) [there seems to be no accepted notation in 
the supergroup case for the distinction analogous to that 
between 0(iV) and SO(A^), nor for that between SO(iV) 
and its covering group Spin(A'^), which we ignore here]. 
Furthermore, the representations R, R are reducible un- 
der S0(2n); they split into two nonisomorphic irreducible 
spinors, each of dimension 2"~^. These two spinors cor- 
respond to even and odd numbers of fermions in the un- 
constrained representation (see App. ^). However, under 
0(2n), i?, R do not split; 0(2n) has irreducible represen- 
tations of dimension 2" [this is related to the fact that 
0(2n) is not a direct product of S0(2n) with Z2, unlike 
the case of 0(2n +1)]. Thus we will still call the models 
spin chains, since they involve irreducible representations 
of their symmetry group, 0(2n). 

When we consider the corresponding nonlinear sigma 
models, via the usual correspondence, we naturally con- 
sider the orbit of the lowest weight in R under 0(2n). 
Due to the disconnected nature of 0(2n), as opposed 
to S0(2n), this orbit 0(2n)/U(n) falls into two discon- 
nected pieces, which are both of the form S0(2n)/U(n) 
as manifolds. Similar statements hold for the superman- 
ifolds in the SUSY formalism. 

In a recent paperEil on the class D of random matrix 
problems, which is the same symmetry class as the RBIM 
fermion problem we are considering, it was emphasized 
that the target manifold of the nonlinear sigma model 
that describes it is 0(2?T,)/U(n) in the replica formalism, 
which has two connected components, corresponding to 
those of the group 0(2n) [or the corresponding super- 
group OSp(2n|2n)]. This opens a possibility not usu- 
ally considered for nonlinear sigma models, that the con- 
figurations include fluctuations (i.e. domains) where the 



sigma model field is on different components of the target 
manifold. This implies that additional parameters, be- 
yond the usual couplings like g,9 ior continuous deforma- 
tions of the field, must appear in the model, to describe 
the domain walls; for example, a fugacity per unit length 
of domain wall. When the fugacity is small, there are 
essentially no domain walls, and the model would reduce 
to that with target space S0(2n)/U(n). 

In our approach, we have arrived at spaces of states 
that correspond to both parts of the target manifold, and 
further the spin chain Hamiltonians contain in general 
eight parameters. Therefore, our spin chain models de- 
scribe a strong-coupling version of the physics of the non- 
linear sigma model with domain walls included. These 
models include the pure Ising model, and weak-disorder, 
limits. Note that the latter are not accessible simply as 
the strong coupling, g — s- 00, limit of the S0(2n)/U(n) 
nonlinear sigma model (compare Ref . |l^) . 

What we have found in this paper is that the states in 
the SUSY description, can be viewed not only as domains 
of two "phases" , but that the discrete (Ising-like) degree 
of freedom, which labels which phase [component of the 
target manifold, or irreducible spinor of S0(2n)] a point 
in 2D space is in, can be replaced by additional continu- 
ous variables. These continuous degrees of freedom turn 
the model into a nonlinear sigma model with symmetry 
S0(2n -I- 1) [or SUSY osp(2n + 1 1 2n)] broken by certain 
terms in the action, or else a strong-coupling version of 
this, at least near the NL. This may be of future use 
in uncovering the physics of these general class D prob- 
lems, not only the RBIM. The replicated spin chains for 
0(2n) at nonzero n have not been considered previously 
(except for the n = 1 case, the usual XXZ model), and 
are also of interest in their own right. Note finally that 
in our earlier discussion of S0(2n -I- 1)- and S0(2n -I- 2)- 
invariant models, the representations of the stated groups 
were irreducible, the corresponding target manifolds were 
connected, and no analogous domain walls were possible. 



IX. CONCLUSION 

In this paper we applied the supersymmetry (SUSY) 
method to analyze an Ising model with a binary distribu- 
tion of random bonds (RBIM). The Nishimori line (NL) 
on the phase diagram of the model is a line with the en- 
hanced SUSY osp(2ri -I- 1 1 2n). On the rest of the phase 
diagram the model has only osp(2n | 2n) SUSY. The en- 
hanced SUSY on the Nishimori line allows us to rederive 



the identities (L4) among various correlation functions. 
More generally, we have shown that the transition on the 
NL has very strong analogies with the integer quantum 
Hall effect transition, and other random fermion prob- 
lems in 2D, such as the spin quantum Hall transition, 
which can also be modeled by (super-)spin chains with 
alternating dual irreducible representations at the sites, 
and staggered couplings. The conformal field theories of 
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the critical points are mostly unknown at present. We 
emphasize that, in view of our results and those of Ref. 
1^, the fixed-point conformal field theory of the multicrit- 
ical point in the RBIM with a generic distribution for 
the bonds (not only those satisfying the Nishimori con- 
dition) must have at least osp(2n + 1 \ 2n) SUSY, and 
this is a requirement for any future proposal for a con- 
formal field theory of the multicritical point within the 
SUSY formulation. We have also demonstrated that such 
higher SUSY points occur in other problems, such as a 
ID model, and probably elsewhere. After analyzing the 
phase diagram of generalized Hamiltonians with the same 
enhanced SUSY as the NL, we suggested that the tran- 
sitions in many or all of these more general 2D mod- 
els are in a universality class with a still larger SUSY, 
osp(2n -I- 2 I 2n). This hypothesis is supported to some 
extent by the weak-coupling RG analysis of the nonlin- 
ear sigma models that correspond to the spin chains. 

Fitting our results into the framework of random ma- 
trix ensembles for such problems is an outstanding chal- 
lenge. It is interesting that the nonlinear-sigma-model 
target manifold we obtain on the NL is (except for a; = 0) 
not in the list of those known to correspond to random 
matrix ensembles in Ref. Possibly there is another 
random matrix theory with special symmetries as on the 
NL. 

There are of course a number of other outstanding 
problems, even for the RBIM. We have hardly touched 
the region below the NL, which remains mysterious. The 
fixed point at K = oo (zero temperature) and p = Pc 
is of particular interest. In this region the system can 
be viewed as a superspin chain, since it is a chain of ir- 
reducible representations of its supergroup, OSp(2n|27i), 
to which the larger SUSY, 0Sp(2n + l\2n), is broken by 
superspin anisotropy terms, similar to the XXZ model. 

Note added: Another numerical work on the multiciit- 
ical point of the ±K RBIM has appeared very recentlyc3. 
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APPENDIX A: EQUALITIES FOR 
CORRELATORS 

We will show in this Appendix that the enhanced su- 
persymmetry present on the Nishimori line in our formu- 
lation allows us to reproduce the results of the type of 
Eq. (O). 



We use the formulation of the correlators in the 
Ising mo de] in terms of paths and modified partition 
functionstSEa. Namely, for a correlator of two spins Si-^ 
and we join the points ii and i2 by an (arbitrary) 
path on the lattice, shift all the coupling constants K 
by i7r/2 along the path, and calculate the modified par- 
tition function Z^™"'^^ for the system with the modified 
couplings. Then the correlator is 



{Sii Si.-; 



^ ^(mod) 



z 



(Al) 



where I is the length of the path. 

In the quantum formalism the vertical coordinate on 
the original square lattice plays the role of imaginary 
time r, and the partition function is given by the super- 
trace of an imaginary time ordered evolution operator 
U, composed of the transfer matrices T^i and T„i for all 
the bonds in the model. Because of the supersymmetry 
the partition function equals 1 by construction (see Eq. 
( p. 26 ) and following) for any realization of the random 
couplings: 



^SUSY- STrT,C/ = l. 



(A2) 



When calculating the correlator (Al) we have to modify 



the couplings in the transfer matrices along the path only 
for one particular replica, say, the first fermionic one. 
Then the correlator will be 



(Si.Si,) = (-z)'STrT,[/("-'i). 



(A3) 



Similarly, when calculating (S'i^S'ij)^™"^, we have to 
modify the couplings for 2m — 1 different replicas, in 
which case we must have 2n > 2m — 1. 

Let us see how the transfer matrices are modified, when 
we shift the couplings by i7r/2. Start with the horizontal 
transfer matrix, assuming that the couplings are modified 
for 2to — 1 fermionic replicas: 



T, 



(mod) 



hi 



exp I 2Ki^i+xXsi + 2^ 2xai I 

2m-l 



(A4) 



Upon averaging over the randomness this becomes 

2 m- 1 



^(mod) _ .2m- 
^2i — * 



(A5) 



Then for the correlator of two spins in the same horizon- 
tal row we have 



2-t-r— 1 2m — 1 



= STrT.U W W (2x„fc), (A6) 

where V ^[U]. 



k—i a—1 
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As before for single transfer matrices, we can rewrite 
the last expression in terms of operators, acting in the 
space T' , using the substitution rules obtained above in 
Sec. En: 



z-\-rx / 



\2m-ll 



z-f-r— 1 ^ 

k=i 



(A7) 



Now comes the crucial point. On the Nishimori line the 
zeroth fermion is supersymmetric with the rest of the 
replicas, so we can replace all XQk in the last expression 
by, say, X2m,k- 



[{SiSi 



i+r— 1 2m 



-\-rx I 



\2rn- 



^1] = STrT.y n 11(2^° 



(A8) 



Then we can safely go back to the original space J- , in 
which the last expression is easily identified as 



[{SiS, 



\27nl 



(A9) 



which proves the relation (1.4) for this particular case. 

Now see how vertical transfer matrices are modified. 
When we modify the coupling for the fermionic replica 1 
on a vertical bond, the vertical transfer matrix for this 
replica is modified from 



cosh 



exp \-2Ki,i+yni, 



(AlO) 



to 



i — ^ exp ( 2Ki,i+yni, 

cosh Ki^i+y 



= 2isinhA'i,i+y exp (2Ki^i+ynii^ . (All) 

(since when K is shifted by i7r/2, the dual coupling K 
changes sign). Adding the rest of fermionic and bosonic 
replicas, we obtain 



ji(mod) 



i(tanh/fi^i+j^)^^»-2"-+i. 



(A12) 



If we modify the coupling for replicas 1 through k, we get 
similarly 



ji(mod) 



{'"{ta.nh Ki^i+y 



,A's.-ELi(2'^ 



-1) 



(A13) 



To average this expression, we have to distinguish the 
cases of odd and even k. For an even k ~ 2m we get 



(mod) 



i2™(tanhif)^s»-E==i(2"--i) 
X (l-p + p(-l)*-) 



2m 



n : 



(A14) 



where we introduced 



. _ 2K'(2r. 



-1) 



The corresponding operator in J^' is 

2m 



rpl{mod) _ ■2r, 



'T[, n 

a=l 



(A15) 



(A16) 



Then for the correlator of two spins in the same column 
we get 

j+r—l 2m 
j+r-1 2m 

= STrT.y W n^-M- (Al^) 

r — J Q — 1 

For odd number k ~ 2m — 1 we obtain instead upon 
averaging 



(mod) _ ■2m-lf^^^-^^ ^^^^,_^2,„-l(2,^^_^) 



T|,;"""^ = 2^'"-^(tanh7^) 

X (l-p + p(-l)*-^'+i 



(A18) 



The last factor here is different from the similar factor in 
Tii. Now it gives 1 when Nsi is odd, and 1 — 2p = e~^^ 
when Nsi is even. In the space J^- this factor may be 



written as e^^ ^\ and the Eq. (A18) is replaced by 

2m- 1 



„/(mod) _ ■2m-lTi/ 



a=l 



where 



„2L*(2noi-l) 



(A19) 



(A20) 



2m- 1 



For the vertical correlator we obtain now 

j+r~l 

[{SiSi+ryf'^'-'] = STvTrV Yl yo.(T) n 

T—j CK—l 

(A21) 

On the Nishimori line due to the enhanced supersymme- 
try (and the fact that L* — K*) the factor j/oi may be 
replaced by y2m,i and we again get the equality of the 
type of Eq. ( |e4| )'. 

The structure appearing in the formulation above for 
the correlators is multiplicative in the bonds which are 
modified along a path connecting the spins. Then it is 
straightforward to generalize the arguments of this ap- 
pendix to the case of arbitrary spin correlators. 
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APPENDIX B: REPRESENTATIONS R AND R 

In this Appendix we review the construction of the 
representations R and R of osp(3|2) in terms of uncon- 
strained fermions and bosons (for details see Ref. 25). 
We also discuss how to form a graded tensor product of 
such representations and obtain the invariant product of 
the superspins G and G. 

To construct the representation R we need only one 
complex boson b and one complex fcrmion / and their 
conjugates 6^, with usual commutation relations. In 
terms of these the generators of osp(3|2) are constructed 
as follows. For an orthonormal basis of the Cartan sub- 
algebra we use 

... J. ,„^±{Pf-l). (Bl) 

In the distinguished system of simple roots of osp(3|2) 
one root ai is odd ( "fermionic" ) , and one root a2 is even 
("bosonic"). The generators corresponding to these roots 
(and their negatives) are 



1 ft 



(B2) 



The other roots are ^ ai + a2, on = ai + 2a2 (both 
odd), a5 = 2ai + 2a2 (even), and their negatives. The 
corresponding generators are 



{-iT^b\ 

b\f\ e_ 
(6t)^ 



bi^iy 



fb, 
--b\ 



(B3) 



Note that the generators corresponding to the roots a2 
and as contain expression (— l)"-f . This is a "twist" op- 
erator for the fermion, which means that it anticommutes 
with / and /t. It is necessary to ensure that these gen- 
erators obey the (anti-)commutation relations. In other 
words, these choices reflect the grading appropriate for 
osp(3|2), instead of that which is natural in the present 
Fock space. 

The vacuum for bosons and fermions |0), defined in the 
usual manner 



6|0) = /|0) = 



(B4) 



is the lowest weight state of the R representation. The 
remaining states are obtained by the action of the raising 
generators, and it is easy to see that they span the whole 
Fock space of b and /. The weights of the states in R in 
terms of and are shown in Fig. 0. We also show 
in this Figure the organization of the states in doublets 
under the gl(l|l) subalgebra generated by 



1 



i? = 716 + 71/, N=-{nb-nf), 
= = 6t/, F = e_„, = fb, 




FIG. 7. The weights of states in R. The action of the 
positive root generators is shown by arrows. The states are 
grouped in pairs which are the doublets under gl(l|l) (see 
text for details). 



(see Ref. ^ which contains detailed discussion of the 
irreducible representations of gl(l|l).) The doublet of 
states with E = m is denoted by Dm- 

From Fig. |^, we can see that the grading of states, con- 
sistent with that of the SUSY generators, and such that 
the vacuum (lowest weight state) is even, is that states 
are even or odd according as the number of bosons is 
even or odd. This agrees with the choice we made in Sec. 
U for other reasons. We may also note that the gener- 
ators without strings, which are bilinears in the bosons 
and fermions, generate the osp(2|2) subalgebra. The lat- 
ter algebra is consistent with the natural grading on the 
Fock space. This is not in contradiction to the above con- 
struction, because the Fock space decomposes into two ir- 
reducible "spinor" representations of osp(2|2), which are 
connected to each other only by the shortest roots e±2, 
e±3 that are not present in osp(2|2). 

The construction of R is similar. The difference is that 
we start with negative norm bosons b and b^ satisfying 



[6,5t] = -l, 



(B6) 



and another pair of the usual fermionic operators / and 
P . One possible choice of the generators of osp(3|2) in 
the R representation is 



hi 



-0-3 



-■0:4 



bp, e_„l=/6^ 

-(-1)"^/, e_„, = -/t(-l)"^ 

6(-l)"N e_„3 = (-l)"^6t. 



-bf, 
12 



-Q4 



= -Pb^, 



Cqs — b , 6— Q5 — (b^) J 
Here the number of b bosons is defined as 



-V^b. 



(B7) 



(B8) 



(B5) 



The minus sign in this expression implies that 71;, 
is a nonnegative integer, with eigenstates \ni) = 



24 



(nb!)~^/^(6''')"5|0). Note that now the 02 andjaa gen- 
erators contain a twist operator for the boson b, which 
also ensures the proper (anti-)coinmutators. 
Now the vacuum |0) for / and b defined as 



where 



/|0> = 5|0> = 



(B9) 



is the highest weight state of the R representation, and 
the remaining states span the whole Fock space of / and 
b. The states of R are now organized in doublets Dm of 
the gl(l|l) generated by 



E = -ni - nj, 



N : 



— (rif 



(BIO) 



Next we have to combine the representations R and 
R in the alternating fashion, as in Fig. H When we try 
to do that we immediately realize that the twist opera- 
tors of individual R and R representations are not ade- 
quate for their job in the tensor product. They should 
be replaced by "strings" , similar to the ones used in the 
Jordan- Wigner transformation. One possible convenient 
choice of these strings is the following. For the represen- 
tation Ri (numbering as in the Fig. |^) the twist operator 
is replaced by 



s.=n(-i) 



(Bll) 



k<i 



k>i 



Similarly, for Ri the operator (— 1)"^; should be replaced 
by 

% = ]^(-i)»^'«+"5<=. (B12) 

k<i k>i 

Note that in fact, — Sj+i- 

For the purposes of Sec. M we need to consider only one 
pair of the antiferromagnetically coupled superspins. In 
this case the E operator common for both representations 
is S = (-1)"/+"'.. First we consider the fully osp(3|2) 
invariant product appearing in Hk on the NL. In terms 
of the root generators this product is given by 

str GG = - hihi - - ^ KZl {eaS-a + Sae-a) , 



a>0 



(B13) 



where K^a — K{e-a, Ga) are the values of the Killing 
form on the pairs e_Q,, Cq. For osp(3|2) they are 

K — —2 K — —4 K — —4 

K^^, = -2, K_^, = -4. (B14) 

With the mentioned expression for E this gives 
- 4 str GG = (6t)2(6t)2 + 52^2 ^ 2{b^b^ + bb){p p + //) 

— 2ribn^ — nil — riij + 2nfnj — Hf — njr 

-iPP + ff + b^P +bb) 



= J^ - J, 



(B15) 



J = PP + If + bm+bb. 



(B16) 



Note that the term J in Eq. (B15) comes from the 
roots ±a2 and ±03, and comes from all the remaining 
roots. These remaining roots are exactly the roots of 
osp(2|2). Therefore, the term is the osp(2|2) invariant 
product. That observation allows us to write the general 
anisotropic product as 

4strAGAG = AJ- (B17) 



APPENDIX C: PERTURB ATIVE BETA 
FUNCTION FOR SO(2ri + 1)/U(n) 

In this Appendix we derive the perturbative beta func- 
tion of the weakly-coupled nonlinear sigma model on 
S0(2n -|- 1)/U(n) target space, to one-loop order. The 
underlying ideas for a general sigma model have been 
discussed extensively by FriedarO, and we can be brief 
(similar calculations can be found in Rcf. We then 
discuss the resulting flows at n = 0. 

Consider a general homogeneous space G/H, where H 
is a subgroup of the group G. The neighborhood of points 
gH S G/H of the "origin" O = eH [g is an arbitrary 
element, and e is the identity, in G) may be parametrized 
in terms of dimG — AimH coordinates by writing 
g = exp{X^Ti} (repeated indices summed). Here Tj 
denotes a basis of the vector space Q /H spanned by the 
generators of G which are not generators of H {Q and H 
denote the Lie algebras). The sigma model on G/H is 
then defined by the action 



S 



d\rnj{X{r)}d^X'{r)d^X-'{r) (CI) 



where r is the coordinate of two-dimensional space. The 
metric rijj{X} on the target space G/H of the sigma 
model serves as the coupling constant (s). Every point 
P = gH in the coset space can be reached from the 
origin by left multiplication, and every element of the 
tangent space at P can be similarly obtained from the 
tangent space G/Ti. at the origin. Therefore, the metric 
at any point P is uniquely determined by that at the ori- 
gin O, where it represents a symmetric bilinear form on 
the vector space G/7i. In order for this bilinear form to 
represent a metric at the origin it must be invariant under 
the subgroup H (which acts by conjugation). Therefore, 
the metrics on the homogeneous space are in 1-1 corre- 
spondence with iJ-invariant symmetric bilinear forms on 

g/n. 

For sigma models on general manifolds (not necessarily 
homogeneous spaces) the two-loop beta function i^ZI 



dl 



1 



Rij{X) + -Ri 



klr 



,{X) i?/'™(X) 



(C2) 
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where X'^ is any system of local coordinates, and 
R''kim{X), Rij{X) are the Riemann and Ricci tensors, 
respectively, at the point of the manifold with coordi- 
nates X . 

For a homogeneous space G/H, it is enough to com- 
pute the beta function for the metric at the origin O = 
eH (since all other points can be reached by left multipli- 
cation with elements of G, acting as isometries), where 
it reads 



dr]ij 



1 



KLM 



(C3) 



This form of the beta function is convenient since it does 
not require reference to any parametrization of the coset 
spaceu-j Rather, the Riemann tensor of the homogenous 
spacecJ'Lj, viewed as a Riemannian space, has a simple 
expression in terms of of the structure constants fij^ 
fjj'^ of the Lie algebra and the metric: 



Rklij 



1 



MM' 



{flML — ImLI + fLIM)^ 
X {fjKM' — IkM'.J + Im'Jk) 
+ -^{fjML — JmLJ + /l./m)?7^^^^ 



IKM' 



Ikm'i + Jm'ik) 



+ 2^1-^ {Jmkl — Jklm + Ilmk) 



(C4) 



where indices K,L,M,M' denote generators in G/Ti-, 
which are lowered and raised by means of the metric rju 
and its inverse rj^"^ . Indices a denote generators in Ti. 
The Ricci tensor is obtained, as usual, by contraction. 



Rl, 



v'^' Rklij- 



(C5) 



We now discuss the space of all possible (G-invariant) 
metrics on the homogeneous space, that is all i?-invariant 
symmetric bilinear forms on the vector space G /TL. Since 
the latter transforms in a (real) representation of H 
(under conjugation), the bilinear form rju — rj{Ti,Tj) 
must, by Schur's lemma, be a multiple of the unit matrix 
on each irreducible component (assuming, for simplicity, 
that each such component occurs only once). Consider, 
for example, the homogeneous spaces SO(A^)/SO(iV — 
1), the familiar 0{N) vector models. Here G /"H = 
so{N)/so{N — 1) transforms in the (irreducible) vector 
representation of SO(A^— 1) and therefore there is only a 
one-parameter family of metrics. This is the case for all 
symmetric spacesEa, whose sigma models have therefore 
only a single coupling constant (the scale of the metric). 

The case of interest in this paper is G/H = S0(2n -|- 
1)/U(n), which is not a symmetric space. It has a 
two-parameter family of metrics, and the corresponding 
sigma model has therefore two coupling constants. To 
see this one notes that the vector space so(2?i -I- l)/u(n) 
decomposes (over the real numbers) under the adjoint 



action of U(n) into two irreducible representations. One 
of them is of dimension n{n— 1); the corresponding gen- 
erators will be denoted T/^ . The other is of dimension 
2n, and the corresponding generators will be denoted Tj^ . 
These vector spaces may be identified with the cosets of 
Lie algebras so(2n)/u(n) and so(2n -I- l)/so(2n), respec- 
tively. This decomposition corresponds to the chain of 
subalgebras, u(n) C so(2n) C so(2n-|-l). The two metric 
components can be specified as follows. Consider first the 
("standard") Cartan-Killing metric K on the entire Lie 
algebra G — so{2n + 1). We choose the basis of genera- 
tors Ti such that K{Ti, Tj) oc 6ij (the structure constants 
with indices lowered by this metric are then totally anti- 
symmetric) . By restriction this is an iJ-invariant bilinear 
form on the subspace G/Ti., on which it is block-diagonal 
on the two irreducible representation spaces of H . The 
scales of the metric on the two blocks represent the two 
parameters of the metric, say 771 > and 772 > 0, and we 
can write explicitly 

VI, J = m^iJi^J,JiK(Ti^,Tj^) +j]2Sij2Sj.j^K(Ti^,Tj^). 

(C6) 

Note that one may relate the structure constants fijK 
of Eq. (C4), with indices lowered with the metric rjjj, 
to those with indices lowered with the Killing metric 
K{Tj,Tj), which are totally antisymmetric. 

The computation of the Ricci (and Riemann) tensor 
of the homogeneous space so( 2n - I- 1 Vu(n) is tedious but 
straightforward, using (C4), ( |C5| ), (C6). In terms of the 
following parametrization of the metric. 



1 

m = 

9 



1 

2xg' 



(C7) 



one obtains from (C3) the one- loop beta functions: 



§ = 2(n-l)5x(.-l)[l-^- ' 
al n — 1 



(C8) 

x]+0{g'). (C9) 



These equations are valid in the limit g ^ with x fixed. 

The parameter a: > measures the relative strength of 
the two metric components. There are two special cases, 
X = and a; = 1, which we now discuss in turn. Consider 
the chain of vector spaces (Lie algebras) u(rt) C so(2n) 
C so(2n+ 1). As a; — > 0, one sees from Eq. (pTl) that the 



stiffness of the fluctuations of the sigma model (CI) asso- 
ciated with the metric component 772, that is of those in 
the space so(2r7 + l)/so(2rt), becomes infinite. At x — 
these fluctuations in the gradients (with respect to r) 
of the sigma model field are forbidden, and the only re- 
maining fluctuations are those associated with the metric 
component 771, that is of those in so(2n)/u(n), together 
with a degree of freedom on SO(2n-|-l)/SO(2n) (a sphere, 
5^") which is independent of r and is therefore global. 
This is related to the structure of SO(27i-|-l)/U(7i), which 
[because of the chain of subgroups U(7i) C SO (277) C 
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S0(2n + 1)] can be viewed as a fiber bundle with base 
space S0(2n+ 1)/S0(2n)^ S'2«^ ^^^^ S0(2n)/U(n). 
Thus for each point on the sphere 5^", there is a copy 
of the space S0(2n)/U(n) in which the field can fluc- 
tuate locally. Because of the global degree of freedom 
on S*^", there is still a global 80(2?! + 1) symmetry. In 
simple terms, the symmetry is spontaneously broken to 
S0(2n); this does not violate the Hohenberg-Merniin- 
Wagner theorem, which applies for integer n > 1, be- 
cause the coupling I/772 = 0. Neglecting the global de- 
gree of freedom, the line a; = now corresponds to the 
S0(2n)/U(n) sigma model. (These remarks explain why 
only non-negative powers of x appear in the perturba- 
tive beta functions.) This lin e is an invariant of the RG 
flow, and the beta function (C8) reduces to that of t 
S0(2n)/U(n) sigma model, which to three-loop orde 



of X of order 1 is reached when I — Iq {Iq is the logarithm 
of the short distance cutoff, the scale at which xo, go 
are defined) is / — /q — ^/i'^go^o), at which g is of order 
502^0- Then we use the two-loop flows at a; — 1, which 
should be suffici ent accuracy, starting from these values. 
Integrating Eq. (8.4), we find finally for the crossover that 
passes close to the two fixed points at x = and x = 1, 
that g becomes of order one when the length scale e' is 



exp[l/{2goXo) + l/{8gixl)], 



(Cll) 



in units of the short distance cutoff, e'" (the numerical 
factors in the exponent should not be taken too literally). 
The two-loop corrections near a; = will generate only a 
factor of a power of g in this length scale. Flows that start 
at a; > 1 and g small give a similar scale, ~ exp[l/(8go)]- 



dg 
dl 



= 2{n - l)g^ + 2{n^ - 3n + 4)5^ 

+ (3n3 - Mn^ + 35?i - 28)/ + O(g^). (CIO) 



For n > 1, g flows to large values. 

At X = 1, on the other han d, o ne can check by direct 
calculation that the metric in (C6) reduces to that of the 
symmetric space SO(2n-|-2) /U(n-l-l) of higher symmetry. 
Therefore, the line x — 1 must also be an invariant of 
the RG flow. On this line, the one-loop beta function 
in (ICSD reduces t o tha t of the symmetric space S0(2n -|- 




Eq. (|C10|) with n 1— > n + 1 , as expected. One 
sees from (C9) that for n > 1 both lines a; = and a; = 1 
are attractive at weak coupling. The line x = (n — l)/n 
is the separatrix between these two regimes. 

There is also a limit, rji — 0, in which fluctuations in 
S0(2n)/U(n) are "soft", and can be gauged away; com- 
pare the discussion in Ref. In this case the model 
reduces to the nonlinear sigma model with target space 
S0(2n+ 1)/S0(2n) = 5*^", mentioned earlier. However, 
this strong-coupling limit cannot be accessed perturba- 
tively in 1/ryi, l/r/2. 

In the replica limit n — s- 0, of interest in this paper, 
we obtain Eqs. (B.S). Note that now, in contrast to the 
case n > 1, a; = is repulsive, and that on the line 
X = 0, g flows towards weak coupling. Near a; = 0, 
the one- loop flow lines are hyperbolas, xg — constant. 
In the vicinity of the line a; = 1 the one-loop flow lines 
are exponentials: g/g* = exp{— 2(a; — 1)}. The one-loop 
flows are qualitatively different depending on the sign of 
{x — 1). When a; < 1 the coupling constant g decreases 
upon RG flow until it approaches some asymptote g = g* 
at X = 1, while for x > 1, on the other hand, g increases 
towards g* and x decreases towards 1, as I 00. 



To two- loop order, as discussed in Sec. VIII , the region 
around x = 1 flows towards strong coupling. We now 
consider the behavior of the flows, in particular those 
which start with bare values near x = 0. Use of the one- 



loop equations (8.3) near x = 0, with bare values xg and 
go, with Xo and goxo assumed small, shows that a value 
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